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Abstract 

Thurston introduced shear deformations (cataclysms) on geodesic 
laminations - deformations including left and right displacements along 
geodesies. For hyperbolic surfaces with cusps, we consider shear defor- 
mations on disjoint unions of ideal geodesies. The length of a balanced 
weighted sum of ideal geodesies is defined and the Weil-Petersson (WP) 
duality of shears and the defined length is established. The Poisson 
bracket of a pair of balanced weight systems on a set of disjoint ideal 
geodesies is given in terms of an elementary 2-form. The symplectic 
geometry of balanced weight systems on ideal geodesies is developed. 
Equality of the Fock shear coordinate algebra and the WP Poisson 
algebra is established. The formula for the WP Riemannian pairing of 
shears is also presented. 

1 Introduction 

As a generalization of the Fenchel-Nielsen twist deformation for a simple 
closed curve, Thurston introduced earthquake deformations for measured 
geodesic laminations. Later in his study of minimal stretch maps, Thurston 
generalized earthquakes to shears (cataclysms), deformations incorporating 
left and right displacements |Thu98] . Bonahon subsequently developed the 
fundamental theory of shear deformations in a sequence of papers [Bon96, 
IBon97a( IBon97b] . At the same time, Penner developed a deformation theory 
of Riemann surfaces with cusps by considering shear deformations on disjoint 
ideal geodesies triangulating a surface [Pe n87} IPen92l IPenl2] . More recently 
shear deformations play a basic role in the Fock-Goncharov work on the 
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quantization of Teichmiiller space [FG071 IFC99[ IFG06j and in the Kahn- 
Markovic work on the Weil-Petersson Ehrenpreis conjecture [KM08 . 

The Weil-Petersson (WP) geometry of Teichmiiller space is recognized 
as corresponding to the hyperbolic geometry of Riemann surfaces. For ex- 
ample, twice the dual in the WP Kahler form of a Fenchel-Nielsen twist de- 
formation is the differential of the associated geodesic-length function. Also 
for example, the WP Riemannian pairing of twist deformations is given by 
a sum of lengths of orthogonal connecting geodesies, see Theorem [3J and 
[Rie05]. An infinitesimal shear on a disjoint union of ideal geodesies is spec- 
ified by weights on the geodesies with vanishing sum of weights for the edges 
entering each cusp. We define the length of a balanced sum of ideal geodesies 
and find that twice the dual in the WP Kahler form of a shear is the differ- 
ential of the defined length. We then present the basic WP symplectic and 
Hamiltonian geometry in Section [7] with Theorem [21] and Corollaries [22| [551 
and [24j The results include new formulas for the Kahler form. We show 
that the Poisson bracket of a pair of weight systems on a common set of 
triangulating ideal geodesies is given in terms of an elementary 2-form com- 
puted from the weights alone. In Section El we use the elementary 2-form to 
show in Theorem [29] that the Fock shear coordinate algebra introduced in 
the quantization of Teichmiiller space is the WP Poisson algebra. The basic 
WP Riemannian geometry of shears is developed in Section [9] with Theorem 
[32l We generalize Riera's WP inner product formula and show that the 
Riemannian pairing of two weight systems on ideal geodesies is given by the 
combination of an invariant of the geometry of ideal geodesies entering a 
cusp and a sum of lengths of orthogonal connecting geodesies. 

There are challenges in calculating shear deformations. In contrast to 
earthquake deformations, shear deformations are in general not limits of 
Fenchel-Nielsen twists and a shear on a single geodesic deforms a complete 
hyperbolic structure to an incomplete structure. For the deformation the- 
ory larger function spaces are involved; for earthquakes geodesic laminations 
carry transverse Borel measures and for shears geodesic laminations carry 
transverse Holder distributions. A general approach would require a de- 
formation theory of incomplete hyperbolic structures. Rather, we follow 
the approach of [Wlp09] and double a surface with cusps across cusps, and 
open cusps to collars to obtain approximating compact surfaces with reflec- 
tion symmetries. Shears are then described as limits of opposing twists. 
Given the above expectations, the approximating formulas include individ- 
ual terms that diverge with the approximation. The object is to show that 
diverging terms cancel and to calculate the remaining contributions. We 
use the Chatauby topology for representations to show that the hyperbolic 
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structures converge and an analysis of holomorphic quadratic differentials 
to show that infinitesimal deformations converge. 

We begin considerations in Section [2] with the variation of cross ratio and 
geodesic-length. A unified treatment is given for Gardiner's geodesic-length 
formula [Gar 75] . Riera's twist Riemannian product formula [Ric05j and the 
original twist-length cosine formula |Wlp83| ■ In Section [3j we review Bona- 
hon's results on shears on compactly supported geodesic laminations and 
Penner's results on shears on ideal geodesies triangulating a surface with 
cusps. The review includes the Thurston-Bonahon Theorem that shears 
on a maximal geodesic lamination are transitive on Teichmiiller space and 
Penner's Theorem on A and h length global coordinate. We include the 
Bonahon-Sozen and Papadopoulos-Penner results that in appropriate set- 
tings the WP Kahler form is a multiple of the Thurston symplectic form. 
In Sections H] and [5j beginning with hyperbolic collars and cusps, we give 
the geometric description of shear deformations and describe the conver- 
gence of opposing twists to shears. In Section [6J we treat the convergence of 
infinitesimal opposing twists to infinitesimal shears. The analysis includes 
the convergence of holomorphic quadratic differentials. In Section UJ we de- 
fine the length of a balanced sum of ideal geodesies and establish the basic 
symplectic geometry results in Theorem 1211 and the following corollaries. In 
Corollary [22l we show that the Poisson bracket of length functions and the 
shear derivative of a length function are given by evaluation of the elemen- 
tary 2-form. We consider the Fock shear coordinate algebra in Section [8l 
We use Penner's topological description of the shear coordinate bracket and 
compute with the elementary 2-form to show that the algebra is the WP 
Poisson algebra. In Section [9] we begin with expansions for gradient pairings 
for geodesies crossing short geodesies. Then in Theorem 1321 we provide the 
formula for the WP Riemannian pairing of balanced sums of ideal geodesies. 
In Example 1341 we calculate the pairing for the Dedekind PSL{2; r L) tessel- 
lation to find an exact distance relation. Finally in Section [10] we give the 
length parameter expansion for the sum of lengths of circuits about a closed 
geodesic. 

It is my pleasure to thank Joergen Andersen, Robert Penner, Adam Ross 
and Dragomir Saric for many helpful conversations and valuable suggestions. 

2 Gradients of geodesic-lengths 

We begin with the basics of deformation theory of Riemann surfaces [Ahl06l 
IHub061 IIT92| . A conformal structure is described by its uniformization. An 
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infinitesimal variation of a conformal structure is described by a variation of 
the identity map for the universal cover. The interesting case for the present 
considerations is for a Riemann surface of finite type, a compact surface with 
a finite number of points removed, covered by the upper half plane HI. For a 
vector field v on the universal cover and parameter e, there is a variation of 
the identity map w e (z) = z + ev + o(e), for z, respectively w, conformal coor- 
dinates for the domain and range universal covers. Provided the vector field 
is deck transformation group invariant, the map is equivariant with respect 
to deck transformation groups. The range conformal structure is described 
by the angle measure aig(dw e ) for the differential dw e = w e ^ z dz + w e ^dz. 
The expansion for the variation provides that dw e = w e ^ z {dz + ev z dz) + o(e), 
and thus arg(dw t ) = arg(u> ej2 ) + avg(dz + ev z dz). The derivative of the vec- 
tor field v z describes the infinitesimal variation of the conformal structure. 
The quantity v z is an example of a Beltrami differential, a tensor of type 

§- z ®Tz. 

For a Riemann surface R of finite type and vector field v defined on the 
surface (equivalently on the universal cover and invariant by deck transfor- 
mations), then w t {z) is a variation of the identity map of the surface and 
in effect describes a relabeling of the points of the surface - the deformation 
is trivial. Nontrivial deformations are given by vector fields on the univer- 
sal cover; vector fields with nontrivial group cocycles relative to the deck 
transformation group. 

We consider B(M), the space of Beltrami differentials on H, bounded 
in L°°. By potential theory considerations, for [i £ B(M) there is a vector 
field v on M with v z = that is actually continuous on H and is bounded 
as 0(\z\ log \z\) at infinity [AB60]. In particular elements of B(M) also de- 
scribe variations of the points of M. We are interested in the corresponding 
variational formula. 

The cross ratio of points of P 1 is given as 

(p-r)(q-s) 

(p,q,r,s) = 

{p-s)[q-r) 

and for q = s + As and rearranging variables, we obtain a holomorphic 
1-form 

(p — q)dz dz dz 

= ( z -p)(z-q) = Jz~-p) ~ Jz^qj- 

The cross ratio and 1-form are invariant by the diagonal action of PSL(2; C) 
on all variables. 

There is a natural pairing of Beltrami differentials with Q(M), the space 
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of integrable holomorphic quadratic differentials on H, 

(p,il>) -> / A# for /i G S(H) and ip G Q(H). 

Rational functions, holomorphic on H, with at least three simple poles on 
R are example elements of Q(H). The holomorphic quadratic differentials 
Q(M) describe cotangents of the deformation space of conformal structures. 
The variational formula for points of R is fundamental. 

Theorem 1. Variation of the cross ratio [Ahl6H [Ahl06| . For p,q,r,s G R 
the variational differential of the cross ratio is 

2 

dlog(p,q,r,s) = — Q pq £l r s G Q(H). 

7T 

The quadratic differentials Q(M) form a pre- inner product space with a 
densely defined Hermitian pairing 

(0,-0)= / (^((fs 2 )" 1 for 0, V> G (5(H) n L 2 



and (is 2 the hyperbolic metric. The pairing is the Weil-Petersson pre-inner 
product [Ahl6lJ WlplO| . The pairing provides formal dual tangent vectors 



for the differentials of cross ratios 

grad log(p, q, r, s) = (dlog(>, q, r, s))(ds 2 ) -1 . 

We are interested for distinct quadruples V = (pi,p%,T\,r%), F = (/i, /2, <7i, 52) 
in the pairing 

(grad log T 5 , grad log F) . 

The pairing is continuous in the quadruples for all points distinct and also 
is continuous for (ri,r2) tending to {pi,P2) and (51,52) tending to (/i,/2)- 
We will evaluate particular configurations for the pairing. 

Let T be the Teichmiiller space of homotopy marked genus 5, n punc- 
tured Riemann surfaces R of negative Euler characteristic. We are inter- 
ested in pairings corresponding to geometric constructions of deformations. 
A point of T is the equivalence class of a pair (R, f) with / a homeomor- 
phism from a reference topological surface F to R. By the Uniformization 
Theorem a conformal structure determines a unique complete compatible 
hyperbolic metric ds 2 for R and a deck transformation group T C PSL(2; R) 
with R = H/r. The Teichmiiller space is a complex manifold with cotan- 
gent space at R represented by Q(R), the space of holomorphic quadratic 
differentials on R with at most simple poles at punctures. 
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The pairing 



(/J,,ip) ->■ / ni> for fi £ B(R) and tp e Q(R) 
Jr 

is the ingredient for Serre duality and consequently the tangent space of T at 
R is B(R)/Q(R) ± |Ahl611 lATuM IHar77l iHubOol HT92 j. The I? Hermitian 
pairing 

= / w{ds 2 r i 

Jr 

is the Weil-Petersson (WP) cometric for Q(R). The metric dual mapping 

<j> -> ^(ds 2 )" 1 G Q(i2) 

is a complex anti linear isomorphism, since Beltrami differentials of the given 
form (harmonic differentials) give a direct summand of Q(R)' L in B(R). The 
metric dual mapping associates a tangent vector to a cotangent vector and 
so defines the WP Kahler metric on the tangent spaces of T\ the mapping 
is the Hermitian metric gradient. 

Geodesic-lengths and Fenchel-Nielsen twist deformations are geometric 
quantities for pairings. Associated to a nontrivial, non peripheral free ho- 
motopy class a on the reference surface F is the length £ a (R) of the unique 
geodesic in the free homotopy class for R. Geodesic-length is given as 
2 cosh l a /2 = tr A for a corresponding to the conjugacy class of A £ T in the 
deck transformation group. Geodesic- lengths are functions on Teichmuller 
space with a direct relationship to WP geometry. A Fenchel-Nielsen twist 
deformation is also associated to a closed simple geodesic. The deformation 
is given by cutting the surface along the geodesic a to form two metric cir- 
cle boundaries, which then are identified by a relative rotation to form a 
new hyperbolic surface. A flow on T is defined by considering the family of 
surfaces {Rt} for which at time t reference points from sides of the original 
geodesic are relatively displaced by t units to the right on the deformed 
surface. The infinitesimal generator the Fenchel-Nielsen vector field t a , the 
differential of the geodesic-length and the gradient of geodesic-length satisfy 
duality relations 

2u>wp( ,ta) = d£ a and equivalently 2t a = Jgrad£ a , (1) 

for ojwp the WP Kahler form and J the complex structure of T (multipli- 
cation by i on B(R)/Q(Q) ± ) |Wlp82j |Wlpi~0l . The factor of 2 adjustment 



to our formulas as detailed in |Wlp07[ §5] is included. 
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We are interested in the WP metric and Lie pairings of the infinitesimal 
deformations grad£ a and t a with geodesic- length functions lp. The formu- 
las begin with Gardiner's calculation of the differential of geodesic-length. 
We now use a single simplified approach that provides Gardiner's dt a for- 
mula [Gar 75] . the cosine formula for t a £p |Wlp83 WlplOj , the sine- length 



formula for t a tp£^ |Wlp83 WlplO| , as well as Riera's length- length formula 



for (grad£ a , grad 1$) [Rie05, WlplOj. The approach combines Theorem [U 



coset decompositions for the uniformization group and calculus calculations. 
An important step is identifying a telescoping sum corresponding to a cyclic 
group action. We present the approach. 

Theorem 2. Gardiner's variational formula [Gar 75] . For a closed geodesic 
a, 

di <* = l E n rAa A (Cz) G Q(R) 
ce(A)\r 

with a corresponding to the conjugacy class of A G T with repelling fixed 
point rA and attracting fixed point a a ■ 

Proof. We begin with the geodesic-length. For a hyperbolic transformation 
A, the geodesic-length is log(As, s, rA, oa) for s a point of R distinct from 
the fixed points. We begin with the variational formula for the cross ratio 
from Theorem [TJ The resulting integrand is in L 1 (H) and H is the disjoint 
union 

|J (J A n C{F) 

nez ce(A)\r 

for T a r fundamental domain. By a change of variables the union over 
domains is replaced by a sum of integrands 



MM = --31 f fiJ2 £ n As s(A n Cz)n rAaA (A n Cz). 
n Jt n ce(A)\r 

The invariance of f2 by the diagonal PSL(2;M) action gives 0, pq (A n w) 
^A~ n pA- n q( w ) an d the given product of forms is 

Q-A-^sA-^siCz) Qr A a A {Cz). 

Using the Q partial fraction expansion, the first factor is 
„ dw dw 



(2) 



A- n + 1 sA~ n s 



(w - A~ n+1 s) (w - A~ n s) 
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and the integer sum telescopes 

N 

^A-^sA-^s = ^A N + 1 sA~ N s 

n=-N 

and as N tends to infinity, A N+1 s tends to oa and A~ N s tends to va- 
(Various forms of the telescoping appear in the calculations for the cosine 
formula |Wlp83j pgs. 220-221], the sine-length formula |Wlp83j pgs. 223- 
224] and the length-length formula [Rie05l pgs. 113-114].) The sum in © 
now becomes the desired sum 

" E V 2 rAaA (Cz). 
c&{A)\r 

□ 

We consider the WP Hermitian pairing of gradients (grad^a, grad^}. 
By ([1]) the imaginary part of the pairing is 

cos 6p. 

pGan/3 

The real part of the pairing (grad^ a , grad^) was first evaluated by Riera 
[Rie05| . We now apply the above approach and with a single simpler treat- 
ment derive the real and imaginary part formulas. Riera's formula involves 
the logarithmic function 



R(u) = ulog 



u + 1 



u-1 



2. 



The function is even with a logarithmic singularity at ±1 and with the 
expansion 

fl( "> = 2 ^ + i + 7^ + -) f0r| " l>L 

In particular for u > 1, the function and its even derivatives are positive 
and the function is 0(u~ 2 ) for u > 1. The function R(u) is also given as 

U 1 Hi 

— tanh -1 — — 2 for Inl > 1 and — tanh _1 u — 2 for \u\ < 1. 
2 u 11 2 11 

We present the pairing formula for the general case of a cofinite group pos- 
sibly with parabolic and elliptic elements. 
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Theorem 3. The complex gradient pairing |Wlp83, Rie05]. For closed 
primitive geodesies a, /3 corresponding to elements A, B £ T, we have for 
the WP pairing 



(gradC grad^) = -5 a/3 e(A)£ a + £ TZ D , 

De{A)\r/(B) 

where 5 a p is the Kronecker delta for the geodesic pair, where e(A) is 2 in 
the special case of the axis of A having order-two elliptic fixed points and is 
1 otherwise, where for the axes axis(A), axis (DBD^ 1 ) disjoint in M, then 

2 

U D = -R(coshd(&xis(A),axis{DBD- 1 ))) 

and for the axes intersecting with angle Qd, then 

2 

TZd = —R(cos6f)) — 2icos9rj. 

Twist-length duality and J an isometry provide that 4(f a , to) = (grad^ a , grades) • 
Proof. For A a hyperbolic element we write 

Ce(A)\r 

and from Gardiner's formula dt a = (2/^)0^ with 

(e A ,e B ) = - f C^lWO ids 2 )- 1 . 

H C£{A)\Y 

We first decompose each left coset (A)\r by considering right (B) cosets and 
then move the (B) action to the two conjugate forms. The resulting sum 
over (B) is telescoping. In particular, we enumerate the cosets of the sum 
by writing for C € (A)\T the decomposition C = DB n , D € (A)\T/(B) for 
n € Z. For A, B primitive hyperbolic elements, we consider uniqueness of 
the presentation of an element of {A)D in the form A m DB n . A non unique 
presentation is equivalent to a solution of A a = DB b D~ 1 for a non trivial 
integer pair (a,b). Since A, B each generate V maximal cyclic subgroups, a 
non trivial solution of A a = DB b D~ 1 provides that A is conjugate to B^ 1 
by the element D. In particular the presentation A m DB n is unique except 
for the case a = (3 with A = DB ±l D~ l . In the case a = /3 we select the 
element A to represent the geodesic and the presentation is unique except 
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for the case of D either the identity or the special case of T containing an 
order-two elliptic E with A = EA~ 1 E. For the special cases there is no 
distinction between left and right (A) cosets; we only use left cosets. The 
special left cosets are for the identity element and the element E. 

Now for each resulting integral of the sum, change variable by writing 
w = B n z; the effect is to move a B~ n action to the variable of VtBss^r B a B - 
Using the diagonal PSL(2;M) invariance of f2, the B~ n action is moved to 
the quadruple of points, resulting in the telescoping sum 

^B n + 1 sB"s^r B a B = ~^r B a B - 

The result is the general formula 

(Oa.Ob) = -s AB ±ie(A) / n Bss n rgaB {ds 2 )- 1 + 

Jm 

y [ n 2 . a . v%~~(ds 2 )-\ (3) 

I r D~ 1 AD a D- 1 AD r BdB y 

De(A)\r/{B) Ja 

where the Kronecker delta indicates that the first integral is only present 
for the case that A = B , e(A) is 2 in the case of order-two elliptic fixed 
points on the axis of A and is otherwise 1, and for the second integral the 
diagonal invariance was used to move the D action to the pair of points. 
For each integral, a change of variable by an element of PSL(2;~R) results 
in the inverse element applied to the tuple of points. It follows that the first 
integral depends only on the PSL(2; R) conjugacy class of B and the second 
integral depends only on the P5L(2;R) class of the pair (D~ 1 AD,B). It 
follows that the first integral is a function of the geodesic-length for B and 
the second integral depends only on the distance between/intersection angle 
of the axes. 

We evaluate the integrals. The differential Q pg is continuous in p, q, 
including at infinity; for q tending to infinity the form limits to dz/{z —p). 
For the first integral of ([3]), we take the pair of points to be and oo, to 
obtain for z = re ie the integral 

f 1 (Bs - s) , . 2 

- / -^l — r, m — 7 r sm 6 rdrdO, 
Ju z z {z - Bs)[z - s) 
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which for P = e Bs, Q = e s becomes 

_ r P {P ~ Q \ - S m 2 9drd9 = 
Jo Jo (r-P)(r-Q) 

(r-Q) 



log 



(r-P) 



■'0 



sin 2 9d9 = -log—. 

2 s 



as expected, since grad£* = 2/7r6*. For the second integral of (|3|), we take 
the first pair of points to be and oo, to obtain the integral 



If ip-q) 



z 2 \ (z — p)(z — q) 
which for P = e ld p, Q = e ie q, becomes 



2 



r 2 sin z 6 rdrd9, 



/•t r°° (p _ n^2 



The r integral has antiderivative 



(P + Q), (r-Q) P Q 

log' 



(P-Q) b (r-P) (r-P) (r-Q)' 

We are evaluating an area integral and 9 varies in the interval (0, 7r); for 
p, q G R, as described, and r real positive, the quotient (r — Q)/(r — P) is 
valued in the complex open lower half plane. The antiderivative is invariant 
under interchanging p, q; we now normalize p to be positive real. We use the 
principal branch of the logarithm; for r close to zero the argument is close 
to —it. Evaluating r at 0, oo and integrating in 9 gives 

— ( — ^— logK — 2) for k the ratio q/p = (q,p, 0, oo). 

To interpret geometrically, compare to [Rie05[ pg. 114], set u = (/c + l)/(re— 
1) = 2(oo, q,p, 0) — 1, to obtain the complex-valued expression 

_(„ log _-2). 

For the lines Ooo and disjoint, the ratio k, = p/q is positive and the 
logarithm is real, with u = cosh 5*, for <5* the distance between the lines. 
For the lines intersecting, the ratio k = q/p is negative and the argument of 
the logarithm is — ir and evaluation gives 

7T 7T 

— i?(cOS#*) — — ITT COS 9^, 
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as desired. 

□ 

The double coset enumeration admits a topological/geometric descrip- 
tion. We consider that a and f3 are primitive and V is torsion-free. On 
the surface R, consider the homotopy classes rel the closed sets a, f3 of arcs 
connecting a to (3. For the universal cover, fix a lifting of a to a line do in 
H; then a connecting homotopy class on R lifts to a homotopy class of arcs 
connecting do to j3 (a line lifting of (5). The relation rel a corresponds to the 
relation of the {A) action on homotopy lifts. In particular, the non trivial 
classes on R rel ot,/3 biject to the classes in EI rel cto,/3, for j3 (disjoint from 
do) ranging over the line liftings of (3 modulo the action of (A); the non 
trivial classes on R correspond to lines f3 disjoint from do- To enumerate 
the pairs (do,/?) for j3 distinct modulo the (A) action, for A generating the 
stabilizer of do and B generating the stabilizer of a line lifting of j3, then line 
pairs distinct modulo the {A) action correspond bijectively to double cosets 
by the rule 

(d ,/3) = (axis(^),axis(Z? J BZ)~ 1 )) corresponds to D € (A)\T/(B). 

The relation axis^-BD -1 ) = D(axis(B)) is part of the correspondence. For 
a finite number of double cosets the corresponding axes intersect. Overall 
the axes enumeration by double cosets, enumerates pairs of line liftings 
of a and j3 modulo the diagonal action of the group T. The geometric 
description comes from the description of a pair of lines. A pair of lines 
either intersects or has a unique perpendicular geodesic, minimizing the 
connecting distance. The cosine and hyperbolic cosine describe the geometry 
of the configurations. 

The present approach to evaluating the pairing is a combination and 
simplification of earlier works. The role of the cyclic group in Gardiner's 
formula was first noted by Hejhal |Hej78[ Theorem 4]. The telescoping 
of the cyclic group sums appears in the proofs of Theorem 3.3 and 3.4 of 
Wlp83| and in Theorem 2 of |Rie05| . although in each case the telescoping 
is presented as a special feature. The basic integral (j4|) is simpler than found 
in the earlier formulations. The present approach can be applied to evaluate 
the second twist Lie derivatives t a tpi^. The first derivative t a £g is a sum of 
cosines of intersection angles. A cosine is given by a cross ratio, the starting 
point for the above considerations. 
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3 Thurston shears 



We are interested in Thurston shears (cataclysms) on ideal geodesies for 
a Riemann surface with cusps. Thurston studied the shear deformation for 
compact geodesic laminations [Thu98]. Bonahon developed the fundamental 
results in a sequence of papers [Bon96, Bo n97aj lBon97bj . We present a brief 
summary of Bonahon's basic results following [Bon96j. In a series of works 
|Pen87l TPen92l lPen!2| . Penner developed a deformation theory of Riemann 
surfaces with cusps by considering shear deformations on ideal geodesies 
triangulating a surface. Our interests include Penner's A-length formulas 
and formulas for the WP Kahler /symplectic form [P P93 ] . We present a brief 
summary of Penner's results following the exposition of the book |Penl2| . 

A geodesic lamination A is a closed union of disjoint simple geodesies. A 
geodesic lamination for a compact surface R is maximal provided R — A is a 
union of ideal triangles. A transverse measure for a geodesic lamination A is 
the assignment for each transverse arc k with endpoints in A c of a positive 
Borel measure fi on the transverse arc with supp(/u) = A n k. If transverse 
arcs kj k are homotopic through arcs with endpoints in A c then the assigned 
measures correspond by the homotopy. The assignment k t— >■ fi(k) is addi- 
tive under countable subdivision of transverse arcs. A measured geodesic 
lamination defines an earthquake deformation by interpreting n(k) as the 
relative left shift of the A complementary regions containing the k endpoints. 
By allowing left and right shifts on complementary regions, Thurston defined 
the shear deformation. The relative left shift of A complementary regions 
again defines a functional on transverse arcs. The functional, called a trans- 
verse cocycle, is only finitely additive under subdivision of transverse arcs. 
A transverse cocycle is not given by integrating a measure, rather is given 
by elements of the dual of Holder continuous functions on transverse arcs. 
The space of transverse cocycles H (A) on a geodesic lamination is a finite 
dimensional vector space. 

Teichmuller space is the space of isotopy classes of hyperbolic metrics. 
A geodesic lamination is represented on each isotopy class of a hyperbolic 
metric. Shear deformations on a given maximal geodesic lamination param- 
eterize Teichmuller space. A projection between leaves is defined for the lift 
of a lamination to the universal covering of the surface. The construction 
begins with the observation that the unit area horoballs in an ideal triangle 
are foliated by horocycles. The tangent field of the partial foliation of ideal 
triangles extends to a Lipschitz vector field on the universal covering; the 
vector field is not defined on the small trilateral regions in each ideal trian- 
gle. The Lipschitz vector field defines a projection between leaves of the lift 
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of the lamination. The projection defines a relative displacement between 
lamination complementary regions. The relative displacement is finitely ad- 
ditive. The relative left displacement is called the shearing cocycle an of the 
surface R. The transverse cocycle for the shear deformation from a surface 
Ri to a surface i?2 is the difference — ur 2 of shearing cocycles. For a 
train track carrying a geodesic lamination, transverse measures are specified 
in terms of non negative weights on the track and transverse cocycles are 
specified in terms of real weights. We also refer to the Thurston symplectic 
intersection form r for a train track. The shearing cocycles for a maximal 
geodesic lamination provide an embedding of Teichmiiller space. 

Theorem 4. [Bon96, Theorems A, B]. The map R i— > ctr defines a real 
analytic homeomorphism from T to an open convex cone C(A) bounded by 
finitely many faces in T~L{\). A transverse cocycle [i is in the cone C(A) if 
and only if r(/x, v) > for every transverse measure v for A. 

The i?-length of the transverse cocycle /i for A is a generalization 

of the total-length of a transverse measure. The i?-length is defined as 



computed locally by first integrating hyperbolic length measure along the 
leaves of A and then integrating the local function on the local space of A 
leaves with respect to the Holder distribution fi. The i?-length generalizes 
the weighted length for weighted simple closed geodesies; i?-length is given 
by the Thurston intersection form and the shearing cocycle as follows. 

Theorem 5. [Bon96, Theorem E]. If fi is a transverse cocycle for the max- 
imal geodesic lamination A and gr £ %(A) is the shearing cocycle of the 
hyperbolic surface R then i^(R) = t(//, o~r). 

The Theorem H] embedding of T into the vector space %(A) provides 
identifications of tangent spaces TT with T-L(X). The identification enables 
a comparison of symplectic forms. 

Theorem 6. [SBOlj . Let R be a compact hyperbolic surface with a maxi- 
mal geodesic lamination A. Then for the tangent space identifications TT — 
T~L{X), the WP Kahler form is a constant multiple of the Thurston intersec- 
tion form. 

A decoration for a hyperbolic metric with cusps is the designation of a 
horocycle at each cusp. Decorated Teichmiiller space T>T is the space of 
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isotopy classes of hyperbolic metrics with cusps and decorations [Penl2j . 
The decorated Teichmuller space is naturally fibered over Teichmiiller space 
with fibers given by varying the horocycle lengths in a decoration. A section 
of the fibration is given by prescribing horocycle lengths. A decoration 
enables a notion of relative length for ideal geodesies. The \-length of an 
ideal geodesic a is A(q) = e s ( a ^ 2 , where 5(a) is the signed distance along a 
between the decoration horocycles; the distance is positive in the case that 
the associated horodiscs are disjoint. We are interested in the A-lengths for 
the isotopy class of a given ideal triangulation A of hyperbolic metrics. An 
ideal triangulation for a genus g surface with n cusps has 6g — 6 + 3n ideal 
geodesies and Ag — 4 + 2n triangles. 




Figure 1: Adjacent ideal triangles with a second diagonal. 

Additional parameters are associated to an ideal triangulation. The ideal 
geodesies divide the decoration horocycles into segments. The h-lengths are 
the lengths of the horocycle segments. For an ideal triangle, the lengths 
are related by h- a = A a /AfeA c , where for the horocycle segment a the triangle 
opposite side is a and the triangle adjacent sides are b, c. We are particularly 
interested in the shear coordinates. An ideal triangle has a median. For a 
pair of triangles adjacent along an ideal geodesic a, drop perpendiculars 
from the medians to a. The shear coordinate for a is the signed distance 
between the median projections; the distance is positive if the projections 
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lie to the right of one another along a. The shear coordinate is given simply 
in terms of A-lengths and h- lengths. In Figure HJ the shear coordinate for 
the diagonal e is given as 



a e = log = log— = log—, [Penl2l Chap. 1, Corollary 4.16]. (5) 



The fibers of the Teichmiiller fibration DT —> T are characterized simply 
by constant shear coordinates. 

By the classical result of Whitehead, triangulations with common ver- 
tices can be related by a sequence of replacing diagonals in quadrilaterals 
[Penl2l Chap. 2, Lemma 1.4]. The effect on A-lengths of replacing diago- 
nals is given by Penner's basic Ptolemy equation A13A24 = A12A34 + A14A23 
for the configuration of Figure [U [Pen 121, Chap. 1, Corollary 4.6]. We also 
note the coupling equation h\hi = h^h^ for the configuration of Figure dj 
the equation follows from the definition of h- lengths. The A and h lengths 
provide global coordinates for T>T- 

Theorem 7. |Penl2t Chap. 2, Theorems 2.5, 2.10; Chap. 4, Theorems 
2.6, 4.2]. For the ideal triangulation A, the X-length mapping T>T — > M> 
is a real- analytic homeomorphism. For V the vertex sectors of the ideal tri- 
angulation, the h-length mapping T>T — > M> * s a real- analytic embedding 
into a real- algebraic quadric variety given by coupling equations. For the 
ideal triangulation, the shear coordinate mapping T — > M A is a real-analytic 
homeomorphism onto the linear subspace given by vanishing of the sum of 
shears around each cusp. The action of the mapping class group MCG is 
described by permutations followed by finite compositions of Ptolemy trans- 
formations. 

The WP Kahler form pulls back to the decorated Teichmiiller space and 
has a universal expression in terms of A and h lengths. We present new 
formulas for the pullback in Section [6j 

Theorem 8. [Penl2t Chap. 2, Theorem 3.1]. For an ideal triangulation A, 
the pullback WP Kahler form on T>T is 



where the sum is over ideal triangles, A* = dlog A* and the individual trian- 
gles have sides a, b and c in clockwise order. 
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The formula is given without Penner's initial 2 factor following the adjust- 
ment to our own formulas as detailed in |Wlp07 §5]. 



Papadopoulos-Penner establish a formula for the pullback ojwp in terms 
of /i-lengths and describe identifications of spaces to establish that 2uj\vp 
coincides with Thurston's intersection form |PP93| . Specifically the au- 
thors show that their change of variable (ff) transforms their formula (f) 
to the formula (fff); the calculation applies to the present setting by tak- 
ing /i(greek index) = — \oghr-r~~ and n(index) = log Xindex and noting the 
factor of 2. 

Corollary 9. [PP93]. For an ideal triangulation A, the pullback WP Kahler 
form is 

ojwp = h a A hp + hp A h y + h y A h a , 
A 

where the sum is over ideal triangles, h* = d log h* and the individual trian- 
gles have vertex sectors a, /3 and 7 in clockwise order. 

In particular the A to h change of coordinates is pre symplectic. 

Papadopoulos and Penner introduce the formal Poincare dual of an ideal 
triangulation. The formal dual is a trivalent graph with an orientation for 
the edges at a vertex. A modification of the trivalent graph is a punctured 
null gon train track. A set of logarithms of A-lengths corresponds to a 
measure on the train track. A modification of the construction of a measured 
foliation from a measured train track parameterizes the space T>A4J- of 
decorated measured foliations. 

Theorem 10. [PP93, Proposition 4.1]. The train track parameterization 
provides a homeomorphism of T>T to T>M.T . The homeomorphism identi- 
fies twice the pullback WP Kahler form and the Thurston intersection form 
2 ojwp = T - 



4 Thurston shears as limits of opposing twists 

We show that weighted Fenchel-Nielsen twists with twist lines orthogonal to 
short geodesies converge to a Thurston shear deformation on ideal geodesies, 
as the short lengths tend to zero. We begin with the collars and cusp 
description [Bus92] . For a closed geodesic a on the surface R of length £ a , 
normalize the universal covering for the corresponding deck transformation 
to be z -> e la z. The collar C{a) = {£ a /2 < argz < vr - £ a /2}/(z -> e la z) 
embeds into R with a the core geodesic. For a cusp, normalize the universal 
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covering for the corresponding deck transformation to be z — > z + 1. The 
cusp region Coo = {Qz > 1/2} /(z — > z + 1) embeds into R. The collars 
about short geodesies and cusp regions are mutually disjoint in R. 

In the universal cover a Fenchel-Nielsen twist deformation for a single 
geodesic line f3 is the piecewise isometry self map of EI with jump disconti- 
nuity across (5 given by a hyperbolic transformation stabilizing /3. A twist 
deformation of magnitude t offsets the /3 half planes by a relative t units to 
the right, as measured when crossing f3. The relative displacement of a com- 
bination of twists on disjoint lines is found as follows. For the displacement 
of q relative to p, consider the twist lines separating p and q (for neither 
point on a twist line). There is a partial ordering of lines based on contain- 
ment of half planes containing p. By definition the (n + l) st line contains 
the preceding n lines in a common half plane with p. The individual twist 
deformations are normalized to fix p. The combined deformation map of EI 
is given by left (post) composition of the individual deformations formed in 
the order of the lines. A basic property is that the Fenchel-Nielsen twists 
on a set of disjoint lines is a commutative group. 

A finite collection of disjoint closed geodesies on a surface R lifts to a 
locally finite collection in H and an equivariant twist mapping is determined 
on relatively compact sets. For our purposes it suffices to analyze finite 
combinations of twists in EI. 

We begin with hyperbolic cylinders and cusp regions. 

Definition 11. For a hyperbolic cylinder with core geodesic 7, an opposing 
twist is a finite combination of weighted Fenchel-Nielsen twists with twist 
lines orthogonal to 7 and vanishing magnitude sum. For a hyperbolic cusp 
region, a Thurston shear is a finite combination of weighted Fenchel-Nielsen 
twists with twist lines asymptotic at the cusp and with vanishing magnitude 
sum. 

A positive shear corresponds to a right earthquake. For a Thurston shear 
an initial piecewise horocycle orthogonal to the twist lines with successive 
displacements given by the negative weights is deformed to a closed horo- 
cycle. The deformed region is complete hyperbolic with a closed horocycle, 
consequently is a cusp region. The vanishing magnitude sum condition is 
required for completeness of the deformed structure. The condition is noted 
in [Bon96, §12.3] and considered in detail in [Penl2, Chap. 2, §4]. 

Lemma 12. The opposing twist deformation of a hyperbolic cylinder is a 
hyperbolic cylinder. The core length of the deformed cylinder is bounded 
uniformly in terms of the initial core length and the twist weights. For 
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a bounded number of bounded weights, the deformed core length is small 
uniformly as the initial core length is small. 

Proof. Opposing twist lines decompose a cylinder into bands, each isometric 
to a region between ultra parallel lines in H. The twist deformation is given 
by translations across lines. The vanishing magnitude sum provides that 
a deformed cylinder is complete hyperbolic containing ultra parallel bands, 
consequently is a hyperbolic cylinder. 

We observe that for disjoint weighted twist lines converging, Fenchel- 
Nielsen twists (normalized with a common fixed region) converge. For a 
core length £, collar twist lines are represented in the band {1 < \z\ < e f ~} in 
H. For £ small, the individual twists are close to the twist line \z\ = 1. The 
magnitude sum vanishing provides that for £ small the combined twist trans- 
formation is close to the identity. In particular for twist weights bounded 
on a compact set the opposing twist is close to the identity uniformly in £. 
The deformed core length is the translation length of z — > ze e conjugated 
by the opposing twists. The deformed core length is uniformly small in £, 
as desired. □ 

Next we make precise the notion of opposing twists converging to a 
Thurston shear and also note a consequence. 

Definition 13. Opposing twists for a sequence of cylinders with core lengths 
tending to zero geometrically converge to a Thurston shear provided the fol- 
lowing. First, the universal coverings are normalized with the hyperbolic 
deck transformations for the cylinders converging in the compact open topol- 
ogy for M to the parabolic deck transformation for the cusp region. Second, 
for a relatively compact open set K in HI whose projection to the cusp region 
contains a loop encircling the cusp, the intersection with K of the weighted 
twist lines for the cylinders converges to the intersection with the weighted 
Thurston shear lines. 



Lemma 14. Consider hyperbolic cylinders converging to a cusp region with 
opposing twists geometrically converging to a Thurston shear. A normaliza- 
tion by HI isometries of the twist deformation maps of HI converges to the 
Thurston shear in the compact open topology for HI. 

Proof. Convergence of lines intersecting a given relatively compact set in HI 
provides convergence on any compact set. As noted convergence of weighted 
lines in HI provides that suitably normalized deformation maps converge in 
the compact open topology. □ 
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Figure 2: A hyperbolic cylinder with geodesies orthogonal to the core 
geodesic and a cusp region with geodesies asymptotic at the cusp. 

5 Chatauby convergence and opening cusps 

The points of Teichmiiller space T are equivalence classes {(R, /)} of Rie- 
mann surfaces with reference homeomorphisms / : F — )■ R from a reference 
surface. The complex of curves C{F) is defined as follows. The vertices of 
C{F) are the free homotopy classes of homotopically nontrivial, non periph- 
eral, simple closed curves on F. A fc-simplex consists of k + 1 homotopy 
classes of mutually disjoint simple closed curves. For surfaces of genus g 
and n punctures, a maximal set of mutually disjoint simple closed curves, 
a partition, has 3g — 3 + n elements. The mapping class group acts on the 
complex C(F). 

The Fenchel-Nielsen coordinates for T are given in terms of geodesic- 
lengths and lengths of auxiliary geodesic segments, [Ab i80 , Bus92, IT92J. 
A partition V = {a±, . . . , a3 9 -3+ n } decomposes the reference surface F into 
2g — 2+n components, each homeomorphic to a sphere with a combination of 
three discs or points removed. A homotopy marked Riemann surface (R, f) 
is likewise decomposed into pants by the geodesies representing the elements 
of V. Each component pants, relative to its hyperbolic metric, has a combi- 
nation of three geodesic boundaries and cusps. For each component pants, 
the shortest geodesic segments connecting boundaries determine designated 
points on each boundary. For each geodesic a in the pants decomposition, 
a twist parameter r a is defined as the displacement along the geodesic be- 
tween designated points, one for each side of the geodesic. For marked 
Riemann surfaces close to an initial reference marked Riemann surface, the 
displacement r a is the distance between the designated points; in general the 
displacement is the analytic continuation (the lifting) of the distance mea- 
surement. For a in V define the Fenchel-Nielsen angle by t? a = 27TT a /£ a . 
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The Fenchel-Nielsen coordinates for Teichmiiller space for the decomposi- 
tion V are (£ ai , , . . . , £ a3g _ 3+n , $a 39 _ 3+n )• The coordinates provide a real 
analytic equivalence of T to (R + x R) 3 3" 3 + n , |Abi801 IBus92l HT92] . 

A partial compactification, the augmented Teichmiiller space T, is intro- 
duced by extending the range of the Fenchel-Nielsen parameters. The added 
points correspond to unions of hyperbolic surfaces with formal pairings of 
cusps. The interpretation of length vanishing is the key ingredient. For an 
£ a equal to zero, the angle i? a is not defined and in place of the geodesic 
for a there appears a pair of cusps; the reference map / is now a homeo- 
morphism of F — a to a union of hyperbolic surfaces (curves parallel to a 
map to loops encircling the cusps). The parameter space for a pair {£ a - ) '& a ) 
will be the identification space R>o X R/{(0,y) ~ (0, y')}. More generally 
for the partition V, a frontier set T(V) is added to the Teichmiiller space 
by extending the Fenchel-Nielsen parameter ranges: for each a G V, extend 
the range of £ a to include the value 0, with $ a not defined for £ a = 0. The 
points of TiV) in general parameterize unions of Riemann surfaces with 
each £ a = 0, a £ V, specifying a pair of cusps. 

We present an alternate description of the frontier points in terms of 
representations of groups and the Chabauty topology. A Riemann surface 
with punctures and hyperbolic metric is uniformized by a cofinite subgroup 
r C PSL(2;M). A puncture corresponds to the T-conjugacy class of a 
maximal parabolic subgroup. In general, a Riemann surface with punctures 
corresponds to the PSL(2; R) conjugacy class of a tuple (r, (Toi), . . . , (Ton)) 
where (roj) are the maximal parabolic classes and a labeling for punctures is 
a labeling for conjugacy classes. A Riemann surface with nodes R' is a finite 
collection of PSL(2;M) conjugacy classes of tuples (r*, (T^), . . . , (r^*)) 
with a formal pairing of certain maximal parabolic classes. The conjugacy 
class of a tuple is called a part of R' . The unpaired maximal parabolic 
classes are the punctures of R' and the genus of R' is defined by the relation 
Total area = 2ir(2g — 2 + n). A cofinite PSX(2;R) injective representation 
of the fundamental group of a surface is topologically allowable provided 
peripheral elements correspond to peripheral elements. A point of the Te- 
ichmiiller space T is given by the PSX(2;R) conjugacy class of a topolog- 
ically allowable injective cofinite representation of the fundamental group 
iri(F) — > r C PSL(2;M). For a simplex a, a point of the corresponding 
frontier space T{a) C T is given by a collection {(T*, (T^), . . . , (Tq^*))} 
of tuples with: a bijection between a and the paired maximal parabolic 
classes; a bijection between the components {Fj} of F — a and the con- 
jugacy classes of parts (T 3 , (T^), . . . , (r^ 3 -)) and the PS'L(2;R) conjugacy 
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classes of topologically allowable isomorphisms TTi(Fj) —> Y J , [Abi771 lBer74j . 
We are interested in geodesic-lengths for a sequence of points of T converg- 
ing to a point of T(c). The convergence of hyperbolic metrics provides 
that for closed curves of F disjoint from <t, geodesic-lengths converge, while 
closed curves with essential a intersections have geodesic-lengths tending to 
infinity, |Ber74l [WIp90l . 

We refer to the Chabauty topology to describe the convergence for the 
PSL(2;M) representations. Chabauty introduced a topology of geometric 
convergence for the space of discrete subgroups of a locally compact group, 
|Cha50| . A neighborhood of T C PSL(2;M) is specified by a neighborhood 
U of the identity in PSL(2;R) and a compact subset K C PSL(2;R). A 
discrete group T' is in the neighborhood Af(Y, U, K) provided Y' n K C YU 
and r n K C Y'U. The sets Af(Y, U, K) provide a neighborhood basis for 
the topology. The PSL(2;M) topology coincides with the induced compact 
open topology for transformations of H. Important for the present consid- 
erations is the following convergence characterization. A sequence of points 
of T converges to a point of T(cr), provided for each component Fj of F — a, 
there exist PSL(2;M) conjugations such that restricted to ni(Fj) the cor- 
responding representations converge element wise to n\{Fj) — > Y J , [Har744 
Thrm. 2]. 

We now consider a Riemann surface R with cusps and data ^ ^>jPj 
for a Thurston shear. The data is a weighted sum of disjoint simple ideal 
geodesies, geodesies with endpoints at infinity in the cusps. The weighted 
sum of segments entering each cusp vanishes. Double the surface across its 
cusps; consider the union of R and its conjugate surface R with the reflection 
symmetry p for the pair. For the geodesic /3j, we write /3j for the union 
Pj U p((3j). To open cusps, given e positive, remove the area e horoball at 
each cusp and glue the remaining surfaces by the map p to obtain a compact 
surface R e . The surface R e has a reflection symmetry (also denoted p) and 
smooth simple closed curves obtained from surgering the j3j (also denoted 
f3 j). The construction provides a homeomorphism from a reference surface 
F to R e for e positive and the simplex a of short curves for F is given by 
the e horocycles. Standard comparison estimates for metrics provide that 
for the uniformization hyperbolic metric, the simplex is realized by short 
geodesies with lengths tending to zero with e. The comparison estimates 
also provide that on the complement of prescribed area collars about the 
short geodesies, the R e hyperbolic metrics converge C°° to the hyperbolic 
metric of R U R, |Wlp90|. The uniformization groups Y(R t ) for the R e , 
Chatauby converge to the uniformization pair Y(R), Y(R), relative to F and 
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the horocycle simplex a. The uniqueness of geodesies and convergence of 
hyperbolic metrics provide that the geodesies /3j in the free homotopy classes 
f3j converge uniformly on a collar complements to f3j U p(/3j) on R(J R. 

We are ready to compare the effect of the Thurston shear ^2bj(f3j U 
—p(/3j)) on RUR to the effect of the opposing twist bjPj on the hyperbolic 
metric of R t . The reflection p reverses orientation and notions of left /right; 
even though R is the mirror image, we require regions to move in the same 
direction by a twist; the minus sign provides the desired effect. Opposing 
twist deformations do not preserve the reflection symmetry. As a preliminary 
matter, we note from Lemma [12] for weights bounded, the opposing twist of 
R e has small geodesic lengths bounded in terms of e. Twisting R e defines a 
family close to the frontier T(cr). We observe the following. 

Lemma 15. For e small and weights bounded, the opposing twist bjfij 
of R e is Chatauby close to the Thurston shear ^ ^jifij U —p(Pj)) of RU R. 
Furthermore, the infinitesimal opposing twist is close to the infinitesimal 
Thurston shear in the sense of infinitesimal variations of PSL(2;M) repre- 
sentations. 

Proof. In brief the convergence of metrics provides for the compact open 
convergence of the twist/shear lines on H, which in turn provides for the 
element wise convergence of representations. By construction of R e , for 
the components Fj of F — a, the representations ir\(Fj) into PSL(2;R) 
converge element wise and the twist lines compact open converge to shear 
lines. Choose generators for the limiting representations and a relatively 
compact open set U C Iff, such that CU D U 7^ for each generator C. 
For e small, the same elements generate the representations of iri(Fj) and 
satisfy the non empty translate intersection condition. The representations 
are completely determined by their action on U. A twist /shear map r 
of EI induces a variation of a representation by varying a transformation 
B by the conjugation tBt^ 1 . Only a finite number of twist/shear lines 
intersect U. The PSL(2;M.) normalized combined twist is given by finite 
ordered compositions as described above. By metric convergence, as e tends 
to zero, on U the twist lines converge uniformly and the twists converge 
uniformly to shears and thus the representations of the finite number of 
generators converge. The representations are element wise uniformly close 
in e. To consider the infinitesimal variations, we introduce a parameter t for 
t bjfij and t ^2 bj(/3jU—p(/3j))- The considerations provide that the initial 
infinitesimal variations of the generators are also close in e. The infinitesimal 
variations of the representations are determined on generators. □ 
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6 Infinitesimal Thurston shears and opposing twists 



We are interested in geodesic-length gradients. A thick-thin decomposition 
of hyperbolic surfaces is determined by a positive constant. The thin subset 
consists of those points with injectivity radius at most the positive constant; 
for a constant at most unity the thin subset is a disjoint union of collars and 
horoballs [Bus92j. Surface representations into P3L(2;M) are Chatauby 
close precisely when their thick subsets are Gromov-Hausdorff close. For 
a sequence of hyperbolic surfaces with certain geodesic-lengths tending to 
zero, we are interested in the magnitude and convergence of geodesic-length 
gradients grad £ a for geodesies a crossing the short geodesic-length collars. 

Applications of convergence of surfaces and gradients include generaliz- 
ing the Gardiner formula, Theorem [21 to balanced sums of ideal geodesies 
and generalizing twist length duality (pQ) to Thurston shears and balanced 
sums of ideal geodesies. The basic matter is to understand the effect of 
Chatauby convergence for sums of the basic differential Q 2 from Section [2j 
We begin with convergence of hyperbolic transformations of HI. 

A hyperbolic transformation with translation length £, fixed points sym- 
metric with respect to the origin and i on its collar boundary is given as 



(i is distance log 1/i to the A axis with endpoints ±l/£). As i tends to zero, 
A converges to the parabolic transformation 



We consider a Chatauby converging sequence of surfaces with short 
length core geodesies and a crossing geodesic intersecting the core geodesies 
orthogonally. A crossing geodesic intersects collars and core geodesies. 
Given a segment of a crossing geodesic a in a thick region, normalize the 
universal coverings so that the segment lifts to a segment along the imag- 
inary axis with highest point at i. Extend the segment by including the 
arcs that connect to core geodesies (the added arcs cross half collars). A 
core geodesic intersecting a lifts to a geodesic orthogonal to the imaginary 
axis. The figures for the universal covers of the surface, Figure HI and the 
Chatauby limit, Figure El are as follows. In the figures the collar lift and its 
limit are shaded. In Figure 01 the left and right circular arcs orthogonal to 
the baseline bound a fundamental domain for a core geodesic transforma- 
tion. 



A 



cosh£/2 l/^sinh£/2 
^sinh£/2 cosh^/2 
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Figure 3: A symmetric compact surface with crossing and core geodesies. 



Chatauby convergence provides that the original segments on the crossing 
geodesic a have length bounded and it is standard that collar boundaries 
converge to horocycles. Figure is the limit of a sequence of Figures H] with 
upper, respectively lower, shaded regions converging to upper, respectively 
lower, shaded regions. The crossing geodesic limits to an ideal geodesic 
connecting cusps. 

Definition 16. For an ideal geodesic a, we write 



for the infinite series, where p,q are endpoints of a lift of a to EI. 

Lemma 17. For a surface R with cusps and an ideal geodesic a, the in- 
finite series d£ a converges. A.s above, consider surfaces R e with reflection 
symmetries obtained by doubling R across its cusps and opening cusps to 
obtain short length core geodesies. Consider that an ideal geodesic a on R is 
approximated on thick subsets by closed core orthogonal geodesies a e on R e . 
There is a Chatauby neighborhood U of RUR such that for R e G U, on thick 
subsets the harmonic Beltrami differentials dl at (ds 2 )~ l and d£ a (ds 2 )^ 1 are 
uniformly bounded and are uniformly close. 

Proof. The d£ a series are bounded by area integrals as follows. We first 
consider regions. In Figured! the unshaded region in H between the shaded 
crescents, by normalization, lies below the line = 1 and outside a circle 
tangent to M at 0. The integral of |S^ooo | 2 = dr/rd9 for z = re t8 over the 
unshaded region is bounded by the integral over the region between the 
shaded sectors in Figure [5] 
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Figure 4: Crossing and core geodesies. The vertical line is the lift of the 
crossing geodesic. The two semi circles orthogonal to the baseline are con- 
secutive lifts of the core geodesic. The left and right circular arcs bound a 
fundamental domain for the hyperbolic transformation stabilizing the larger 
semi circle. The shaded sectors are lifts of half collars for the core geodesic. 
The region bounded by the shaded sectors and the circular arcs covers a 
region containing a component of the thick subset of the surface. 



On a thick region of a surface a holomorphic quadratic differential satisfies 
a mean value estimate in terms of the integral over a hyperbolic metric ball 
of a radius ro at most the injectivity radius. The thick regions of R t and R 
are contained in the projection of the indicated unshaded regions in Figures 
H] and El By the standard unfolding, the absolute values of di ae (ds 2 ) -1 and 
d£ a (ds 2 )~ 1 at a thick point are bounded by the integral of |Oooo| 2 over the 
disjoint union of ro balls about the orbit of the lifted point in the unshaded 
region |WlplO, Chapter 8]. By the above considerations, the integrals are 



uniformly bounded, establishing the first result. 

For the second conclusion, given 5 positive, choose a relatively compact 
set K in the Figure U] region between shaded crescents, such that the integral 
of |^ooo| 2 over the complement between the shaded crescents is bounded by 
5. The sum of evaluations of figoo a * points not in K is bounded by 5 by 
a mean value estimate. Chatauby convergence provides convergence for the 
sum of evaluations of f^ooo f° r the ° rrj ft points in K. Boundedness and 
convergence are established. □ 

Example 18. The ideal geodesic series dt a for a hyperbolic cusp. For a 

cusp uniformized at infinity with integer translation group then the sum over 
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Figure 5: An ideal geodesic and horoballs. The central vertical line is the lift 
of the ideal geodesic connecting cusps. The left and right vertical lines bound 
a fundamental domain for the parabolic transformation stabilizing infinity. 
The shaded sectors are horoballs about the cusps. The region bounded by 
the shaded sectors and the vertical lines covers a region containing the thick 
subset of the surface. 

the group is 



The formula for the integer sum gives d£ a — 27r csc^ ttz dz^ . From the above 
lemma, for a hyperbolic cylinder the series d£ a approximates the cosecant 
squared in the compact open topology o/H. 

We now combine considerations to obtain a uniform major ant for an 
opposing sum of twists and gradients of geodesic-length functions. The 
majorant is the necessary ingredient for general limiting arguments. We 
codify the situation as follows. 

Definition 19. A crossing configuration is a compact surface with reflection 
symmetry with fixed locus a finite union of small length core geodesies 7 and 
no other geodesies having small length. A crossing geodesic a is symmetric 
with respect to the reflection with two intersections with the core geodesies. 
For a crossing configuration, a sum a j^aj of crossing geodesies length 
functions is balanced provided for each core geodesic 7, the weighted inter- 
section number ^ a.j#(aj II7) vanishes. For a surface with cusps, a formal 
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sum a j^aj of ideal geodesies length functions is balanced provided at each 
cusp the weighted intersection number ^ ajf^(ajDh) with each small closed 
horocycle h vanishes. 

Balanced is the precedent to the condition of the weight sum vanishing 
for each cusp for a Thurston shear. To prepare for a convergence argument, 
we first consider the distribution of mass of a harmonic Beltrami differential. 

Lemma 20. A balanced sum a = ^ Oj grad^ Qj . of gradients for a crossing 
configuration is bounded as follows. On the thick subset the absolute value 
\o~\ is uniformly bounded. On a core geodesic 7 collar, uniformized as 1 < 
\z\ < e^ 1 , < 9 < 7r — £j for z = re 1 ® G M, the balanced sum a is bounded 
as 



0((4 + e- 27r ^ + e 2 ^-^)l 



7 

The bounding constants depend only on the number of crossing geodesies, the 
norm of the weights and a choice of Chatauby neighborhood for the limiting 
cusped surface. 

Proof. A general bound for a harmonic Beltrami differential on a 7 collar is 
|/x| is 0((|(^,gradlog£ 7 )| + (e" 2 ^ + e 2 ^ e -^)£- 2 ) sin 2 9 M) (6) 



for M the maximum of \x on the collar boundary |Wlpl2 Prop. 6]. We use 



Theorem [3] to bound the pairings (grad-^ . ,grad£y). By setup the crossing 
and core geodesies are orthogonal. Each core geodesic intersection con- 
tributes —2 to the pairing evaluation. From the balanced hypothesis, the 
weighted sum of intersection contributions vanishes. Each remaining term 
of the evaluation involves a connecting geodesic segment that crosses the 7 
half collar; the width of the half collar is — log£ 7 . For large distance, the 
formula summand 1Z is approximately e~ 2d ^' a \ In [WlplO, Chap. 8] we 



showed that the sum of distances from a to the 7 collar boundary is uni- 
formly bounded. It follows that the contribution I 2 of the half collar width 
can be factored out of each summand. The sum evaluation is 0(£ 2 ), the 
desired bound. Lemma [17] provides the desired bound for a on the thick 
subset. □ 



7 The symplectic geometry of lengths 

There is a length interpretation for a balanced sum A = Yl a j^aj of ideal 
geodesies length functions as follows. Let T~L be a neighborhood of the cusps 
given as a union of small horoballs, one at each cusp. The length L(A) of the 
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balanced sum is the sum with weights Oj of lengths of segments 0£jP\(R—l-L). 
The balanced condition provides that the length does not depend on the 
choice of horoball neighborhood %. For a crossing configuration the length 
of a balanced sum L(A) is defined in the corresponding manner. In the 
crossing case, the value L coincides with the sum of geodesic lengths. 

The length L(A) of a balanced sum is a generalization of the i?-length 
of a transverse cocycle. The balanced condition at cusps is discussed in 
|Bon96[ §12.3], where it is noted that the condition provides a well-defined 
notion of length. The definition in terms of horoballs shows that the length 
L(A) is given as ^ 2oj log X a . for the A-lengths of the ideal geodesies and 
a decoration. An example of a balanced sum is a shear coordinate <7*, see 
formula ([5]) ; the sum is balanced at each vertex of the quadrilateral of Figure 
[H A second example comes directly from the shear coordinates of Riemann 
surfaces. By Theorem[7J the sum ^ crj£ a is balanced since the sum of shear 
coordinates around each cusp vanishes. The adjustment of a factor of 2 to 



our formulas as detailed in [ Wlp07 §5] is included in the following. 



Theorem 21. For a surface R with cusps and a balanced sum A = aj£ aj 
of ideal geodesies length functions, the length L(A) is a differentiable func- 
tion on the Teichmuller space of R with 

dL(A) = J2 a 3 di <*i G 
The formal sum ^ &jCtj is data for an infinitesimal Thurston shear a a with 

The WP twist-length duality 

2oj W p( ,o~a) = dL(A) 

is satisfied. In particular, the Thurston infinitesimal shear o~a is a WP 
symplectic vector field with Hamiltonian potential function L(A)/2. 

Proof. We first observe that L is a differentiable function on the PSL(2; R) 
representation space. For the reference surface F, a simple loop 5 £ tti(F) 
about the cusp has representation into PSL(2; R) a parabolic element that 
generates a maximal parabolic subgroup. Prescribing an area value (at 
most unity) for the quotient of a horoball by the maximal parabolic sub- 
group determines a horoball and horocycle. (The prescription is equivalent 
to a choice of decoration in the Penner approach [Pen04l IPenl2] .) For a 
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pair of elements of tti(F) denning distinct maximal parabolic subgroups, 
the distance between the prescribed horocycles is a smooth function of the 
PSL(2;M) representation. The length L is a sum of distances between horo- 
cycles and hence a smooth function. The differential dL is an element of 
Q(R). In particular the integral of the element over small neighborhoods of 
the cusps is small. The construction of the function and its differential is 
also valid for the distance between collar boundaries. 

Consider a sequence of compact surfaces R t with reflection symmetries 
obtained by doubling and opening the cusps of R. From LemmafTTl on thick 
subsets, the differentials of geodesic-lengths converge uniformly to differen- 
tials for ideal geodesies. From Lemma [20| for a balanced sum, the sum of 
differentials is uniformly bounded in each core collar; the integral of the sum 
is uniformly small over small area collars. As R 6 limits to R, the distance 
between collar boundaries limits to the distance between horocycles. And 
for closed geodesies j3 contained in the thick subsets, the Fenchel-Nielsen 
twists on P U p((3) of R e converge to the twist of RU R and the twist deriva- 
tives of distance converge. The considerations of Chatauby convergence and 
Lemmas [T71 and 1201 can be applied for the Fenchel-Nielsen twists on /3U/?(/3). 
The conclusion is again that the gradient pairing integrals over small area 
collars and small area horoballs are uniformly small. It follows that the 
pairing for a balanced sum length differential and twist converges to the 
limiting pairing as e tends to zero. The derivative of length converges to the 
derivative of length. Reflection-even twists span the reflection-even tangent 
space. The dL formula is established. 

The considerations for infinitesimal Thurston shears are analogous. The 
deformation is smooth and by Lemma [151 the infinitesimal deformation is a 
limit of opposing twists. The opposing twists satisfy ^ a,jt aj = i/2 ^ <Xj grad^, 
on the side of R e that limits to R. We find the e tending to zero limit by 
Lemmas [T7] and [20l The conclusions follow. □ 

We remark that symmetry is basic to considering the R e to R limit of the 
tangent-cotangent pairing. With respect to the reflection p, the differential 
of the length L(A) is even, while the opposing twist and its limit are odd. 
Also the Kahler form is odd since the reflection reverses orientation for 
surface integration. The above duality relation 2ujy/p{ ,&a) = dL(A) is 
established for reflection even tangents of R U R and cannot be applied to 
evaluate a shear pairing u}wp(&b, &a)- 

To evaluate the pairing of Thurston shears, we introduce an elementary 
alternating 2- form for coefficients summing to zero. For a balanced sequence 
{ a j}j=D we consider the partial sums Aq = 0,Ak = Ylj=i a ji 1 < k < p, 
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where by hypothesis A p = 0. We introduce a pairing for balanced sequences 
u(M,{bj}) = ^{Aj + Aj^bj. (7) 

3=1 

We explain that the pairing depends only on the joint cyclic ordering of 
the sequences and that the pairing is alternating. A cyclic shift in the 
index j, 1 < j < p, has the effect of adding a constant to the partial sums 
Aj,0 < j < p. The balanced condition for the sequence {bj} provides that 
the pairing is unchanged. For the alternating property, we have summation 
by parts for balanced sequences {fj} and {gj} with partial sums Fk and Gk 

n—l n 

F k9k+i = F n G n - ^ G k fk- 

k=m k=m 

In particular we have that 

p p— l p 

y Ajbj = ApB p - ^Bjcij+x = - y Bj^aj 
j=i j=i j=i 

and 

p p p 

y Aj-ibj = A p B p - ^Bjcij = - ^ Bjdj 
j=i j=i j=i 

using that Aq,A p ,Bq and B p vanish. The pairing can be written in the 
alternating form 

1 P 

w({Oj}, {&,}) = - A 3 h J - B 3 a 3- (8) 

i=i 

We note that balanced sequences have an interpretation as tangents to the 
regular (p — l)-simplex and oj an interpretation as a closed 2-form on the 
regular simplex. 

The form uj can be evaluated for a pair of balanced sums for a common 
set of disjoint ideal geodesies limiting to a cusp. For balanced sums A = 
^2aj£ a .,B = ^2bj£ a - and a given cusp, consider the geodesic segments 
limiting to the cusp; some geodesies ctj may not limit to the given cusp and 
some may have both ends limiting to the cusp. Choose and label a limiting 
geodesic as the first and enumerate limiting geodesies in the counterclockwise 
order about the cusp. Evaluate the form uj on the enumerated sequences of 
weights {cij} and {bj}. 
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Figure 6: Shear lines at a cusp. The longer horizontal lines represent ideal 
geodesies ending at a cusp on the far left; the uppermost and lowermost 
horizontal lines are identified. The dotted vertical represents a closed horo- 
cycle in the undeformed hyperbolic structure and the shorter solid verti- 
cals form a closed horocycle after applying a shear ajs, for the horizontal 
lines. The shorter verticals are successively displaced by horizontal incre- 
ments — a\, —a>2, . . . , —a p . The shaded horizontals indicate segments along 
the upper and lower edges of each ideal geodesic, segments connecting the 
horocycles of the deformed structure. 

Corollary 22. For the balanced sums A = aj£ aj and B = ^ bj£ a . for a 
common set of disjoint ideal geodesies, the shear pairing is 

u W p(va,o-b) = -o-aL{B) = ^ ^ w({a i },{6 j }), 

~~ cusps 

The Poisson bracket for the length functions L(A) and L(B) is 
{L(A),L(B)} = 2 5>({ aj }, {&,•}). 

cusps 

Proof. The shear-length duality comes from Theorem [2TJ The first line of 
equations is established by finding the contribution to the change in the 
length L(B) from the change in the determination of a closed horocycle at a 
cusp. We refer to the schematic Figure [6] for the basic geometry. To evaluate 
the change in length and uj, geodesic segments are labeled as described 
above. In the a a deformed hyperbolic structure, the distance between closed 
horocycles measured on the upper edge of an ideal geodesic agrees with 
the distance measured on the lower edge. We can compute the change in 
distance by averaging the change for the upper and lower edges. In Figure 
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[6l the change in the first distance is Ai/2, while the change in the j 
distance is (Aj + Aj_i)/2. For the weighted length L(B), the weight for the 
j th distance is bj. The change in weighted distance for the given cusp is 
XX A? + Aj-i)bj/2, as desired. 

We next consider the Poisson bracket. The non degenerate Kahler form 
u>WP defines an isomorphism from tangent to cotangent spaces and a dual 
form uJwp. For the Hamiltonian length functions the Poisson bracket is 
defined as Cowp(dL(A), dL{B)). By duality the pairing is 4luwp(&Ai a s)- 
The final formula follows. □ 

There is a counterpart to Theorem [5] for the setting of shear coordi- 
nates!^] First given an ideal triangulation A, Theorem [7] provides a bijection 
between balanced sum shears X^ a j 5 j an d T as follows, for Sj denoting the 
shear deformations on the A edges. A basepoint R/± £ T in Teichmiiller 
space is determined by all shear coordinates vanishing. The surface R& is 
constructed by gluing ideal triangles with medians on sides always match- 
ing. Each marked Riemann surface R £ T is given uniquely as a balanced 
sum shear or = J2 a j(R) 5 j °f the surface Ra- We show the balanced sum 
length functions are linear in the shear coordinates as follows. 

Corollary 23. For a balanced sum B = X} bj£ a . of lengths of ideal geodesies 
of the triangulation A and a marked Riemann surface R G T then 



Proof. First we observe that all balanced sum length functions vanish at i?A • 
Given a balanced sum B = bj£ aj , consider the double sum of weights 



where the index m enumerates cusps and the index n enumerates half edges 
entering a cusp. The balanced sum condition is the vanishing of the inner 
sums. Each triangulation edge enters two cusps; the enumeration includes 
each triangulation edge twice. Thus the sum of weights of a balanced sum 
vanishes. Since the shear coordinates of i?A vanish, we can introduce a 
decoration T~L for R& such that all /i-lengths have a common value. It follows 

1 Theorem [S] is formulated for left twists/shears while the present results are formulated 
for right twists/shears. The orientation difference explains the interchange of entries when 
comparing 2-forms. 



L(B)(R) = 2>({o,-(ll)},{6,}). 



cusps 






cusps edges at cusp 
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that all A-lengths have a common value Ao- The length L(B) = bj2log Xq 
of the balanced sum vanishes at 

Given a surface R, the path of shears at = t^2 a j(R) s j connects the 
surfaces i?A and R. Corollary [22] provides that the i-derivative of L{B) 
along the path has the constant value ^2 cusps ^>({aj(R)},{bj}). Integration 
in t provides the desired formula. □ 

By Theorem the shear coordinates for the edges of an ideal triangu- 
lation provide a continuous immersion into Euclidean space. In particular 
the shear coordinates for appropriate subsets of edges provide continuous 
coordinates for Teichmiiller space. A procedure determining appropriate 
subsets of edges is given in the proof of Lemma [26] below. From Theorem 
for a subset of shear coordinates without linear relations, the differentials 
of the coordinates are generically linearly independent. Furthermore from 
Corollary [22] for a subset of shear coordinates without linear relations there 
are sets of balanced sum length functions with constant full rank Poisson 
bracket pairing. It follows from the pointwise full rank pairing that the 
differentials of the shear coordinates in the subset are pointwise linearly in- 
dependent on Teichmiiller space. It also follows that the shear coordinates 
from the subset give a basis for the vector space of balanced sums of length 
functions. 

In |Wlp83[ §4], we found for surface fundamental group representations 
into PSL(2; R) that the Poisson bracket of trace functions is a sum of trace 
functions. The present result describes a simpler structure. By construc- 
tion Thurston shears on a common set of ideal geodesies commute and ac- 
cordingly the Poisson bracket of Hamiltonian potential length functions is 
constant. 

We now express the 2-form u in terms of /i-lengths and use the formula 
to give the relation to Corollary [9] 

Corollary 24. For an ideal triangulation A, the pullback WP Kdhler form 
is 

v 

uwp = ^2 hj A hj+i, 

cusps 3=1 

where the first sum is over cusps, the second sum is over h-lengths at a cusp 
enumerated in counterclockwise cyclic order and h* = dlogh*. For an ideal 
triangulation A, the pullback WP Kdhler form is also given as 

ujwp = - ^ d log A e A da e . 
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Proof. We begin with shear coordinates for T and the shear pairing uJwpi^A^ °"b) 
of Corollary 1221 above. The coefficients {a,j}, {bj} are the evaluations of the 
differentials {do~ e } of the shear coordinates on the Thurston shears 04, ag. 
From ([5]) and Figure [TJ the differential of a shear coordinate is d\ogh"/h' 
where h" is the /i-length clockwise from the edge and h! is the /i-length coun- 
terclockwise from the edge. We now write the sum ([7]) at a cusp in terms 
of increments of /i-lengths. We use the notation of formula ([7]). Introduce 
a decoration for the surface and write the shear coordinate increments in 
terms of h- length increments as a,j = hj-\ — hj and bj = gj-i — gj, where 
/i* , g* are now the evaluations of the differential d log h* . The partial sums 
are Aq = and A k = Y2j=i a j = h v ~ hji where with the cyclic ordering 
ho = hp and by hypothesis Y^!j=i hj = 0. We find the contribution to to from 
an individual increment g~k by considering 

(A k + A k ^)b k + (A k+1 + A k )b k+1 = 

(2h p - h k - /i fc _i)(<7 fc __i - g k ) + (2h p - h k+ i - h k )(g k - g k+ i). 

The overall contribution is 1 — h k+ i)g k . We now have that 



and the first formula is established. 

The second formula follows from Theorem [8] and formal considerations. 
From formula ([5]) we have that 



where the ordered side pairs (a, e), (e, b), (c, e) and (e, d) are in counterclock- 
wise order relative to their containing triangles. The pairs are the side pairs 
of Figure [T] with one side a diagonal. Now given a pair of adjacent sides of 
the triangulation A, the pair occurs in two quadrilaterals with one of the 
sides being a diagonal. It follows that the sum of dlogA e A da e over edges 
is twice the sum of Theorem [HJ The second formula follows. □ 

An observation of Joergen Andersen provides a direct relation of the 
above to Corollary [9) The coupling equation /11/12 = /13/14 gives the 2-form 




k=l 




d log A e A da e = A a A A e + A e A Af, + A c A A e + A e A Arf 
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equation hi A h% + ]}2 A h% + /13 A ^4_+ ^4 A fti = for h* = <ilog h*. The 
relation h\ A /t2 + ^3 A I14 = /13 A h<i + hi A /14 follows. Beginning with 
Corollary [9] and referring to Figure [H we observe the following. For an edge 
e of the triangulation, the wedge of h- lengths adjacent to e of the triangles 
adjacent to e can be replaced with the wedge of h- lengths for consecutive 
vertex sectors at the cusps at the ends of e. The replacement agrees with 
the orientations of the formulas. The replacement for each edge of the 
triangulation transforms the first adjacent by side formula to the second 
adjacent by vertex formula. 

Example 25. The form oj for a once punctured torus. 
A choice of three disjoint ideal geodesies decomposes a once punctured torus 
into two ideal triangles. The torus is described by edge identifying two ideal 
triangles to form a topological rectangle with diagonal 7, and then separately 
identifying the horizontal edges a and vertical edges /3. The pattern of 
geodesies at the cusp is twofold a, 7, (5. Consider the triples of balanced 
weights {a, b, —a — b} and {c, d, — c — d} for the sequence a, j3 and 7. For the 
geodesies enumerated according to the pattern at the cusp, the sequence of 
partial sums for the second set of weights is Aq = 0, A\ = c, A2 = —d and 
A3 = 0. The sum (|7|) evaluates to (ca + (c — d)(—a — b) H — db) = (ad — be). 

We now follow the discussion of Bonahon |Bon97bt Theorem 15] and 
Harer-Penner [PH92, Section 2.1] for the dimension of the space of balanced 
sum coefficients. 

Lemma 26. For a surface with cusps and a maximal configuration of dis- 
joint ideal geodesies, the space of balanced sum coefficients has the same 
dimension as the Teichmiiller space. 

Proof. Consider a configuration of ideal geodesies with weights as a graph 
with weighted edges. The graph is connected since ideal triangles fill in the 
configuration to form a connected surface. We will sequentially coalesce and 
remove edges, each time decreasing the number of vertices, to finally obtain 
a single vertex graph. For a surface with a single cusp no coalescing of edges 
is necessary. Otherwise by connectedness, there is an ideal triangle with 
not all vertices at the same cusp. Begin with such a designated triangle. 
If only two vertices are at distinct cusps, then we begin by coalescing an 
edge connecting the distinct vertices. If all vertices are at distinct cusps 
then we begin by sequentially coalescing two edges of the triangle and the 
third edge will not be subsequently coalesced. We label the ends of edges 
as incoming or outgoing at coalesced vertices as follows. Label the ends of 
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edges adjoining the first vertex as incoming. Coalesce the first designated 
edge, remove the weight and label the remaining ends of edges at the second 
vertex as outgoing for the coalesced vertex. At the coalesced vertex the 
weight condition is that the sum of incoming weights equals the sum of 
outgoing weights. To continue, take a path of edges to an uncoalesced 
vertex and coalesce the first edge to an uncoalesced vertex along the path. 
Label the new ends of edges at the coalesced vertex as the opposite type as 
for the initial segment of the coalesced edge. At the coalesced vertex the 
weight condition continues to be that the sum of incoming weights equals 
the sum of outgoing weights. Continue coalescing edges until only a single 
vertex remains. For a surface of genus g with n cusps, there are 6g — 6 + 3n 
edges in a maximal configuration. A total of n — 1 edges are coalesced and 
6g — 6 + 2n + 1 edges remain. At least one edge of the initial designated 
triangle gives rise to an incoming-incoming edge of the final coalesced vertex. 
The single weight sum relation is a non trivial condition for the weight on 
the incoming-incoming edge. The space of weights on the final graph has 
the expected dimension. □ 



8 The Fock shear coordinate algebra 

Fock and Goncharov in their quantization of Teichmuller space introduced 
and worked with a Poisson algebra for the shear coordinate functions [FG071 
FC99l lFG06| . The quantization considerations begin with the Fock-Thurston 
Theorem that for any ideal triangulation, the corresponding shear coordi- 
nates (without the vanishing sums about cusps condition) provide a real- 
analytic homeomorphism of the holed Teichmuller space to Euclidean space 
[Penl2t Chap. 4, Theorem 4.4]. Fock proposed a Poisson structure by intro- 
ducing a natural bivector, an exterior contravariant 2-tensor rj and defining 
{ft 9} = i(df,dg),r]) for f,g smooth functions. A relationship to the WP 
Kahler form was also proposed. A bivector defines a Poisson structure with 
Jacobi identity provided its Schouten-Nijenhuis tensor vanishes. 

Theorem 27. [FG07llFC99| For an ideal triangulation A and corresponding 
shear coordinates, the bivector 

v - d d d d d d 



r]A = > t; — A — h 7: — A — h — — A 



do a da b da b da c da c da a 

is natural for the holed Teichmuller space, where the individual triangles 
have sides a, b and c in counterclockwise order. 
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Penner gave a topological description of the bracket of shear coordinates 
|Penl2l pg. 81], a proof that the bivector is independent of triangulation 
and also determined the center of the algebra |Penl2[ Chap. 2]. For the 
topological description of the bracket, recall the definition of the fat graph 
dual to an ideal triangulation. To construct the fat graph G embedded in 
the surface, choose a vertex interior to each triangle and connect vertices by 
an edge when triangles are adjacent. The result is a trivalent graph with a 
cyclic ordering of edges at each vertex. The trivalent graph is a deformation 
retract of the surface. 

Penner 's topological description of the bracket is the following [Penl2, 
pg. 81]. Consider an ideal triangulation A with dual fat graph spine G. If 
a, b G A are distinct edges, then let e a b be the number of components of the 
complement of AUG whose frontier contains points of a and b, counted with 
a positive sign if a and b are consecutive in the counterclockwise order in the 
corresponding region, and with a negative sign if a and b are consecutive in 
the clockwise order H Setting e aa = for each a G A, e a b takes the possible 
values 0, ±1, ±2 and comprises a skew-symmetric matrix indexed by A. The 
quantity e a b is the count of oriented vertex sectors jointly bounded by a and 
b. 

Definition 28. The Fock shear coordinate algebra is defined by the bracket 
{o- a ,o-b} = tab for a,b £ A. 

From formula ([5]) and Figure [TJ a shear coordinate is a balanced sum of 
length functions. For Riemann surfaces with cusps the WP Poisson bracket 
of sums of length functions is given in Corollary [22] in terms of weights and 
the form uj. We evaluate uj for quadrilaterals and find that the evaluation 
agrees with Penner's topological description of the count e a b- 

Theorem 29. The Fock shear coordinate algebra is the WP Poisson algebra. 
The Fock shear coordinate bracket is given by the form uj. 

Proof. We begin with Corollary [22] providing that the Poisson bracket of 
the shear coordinates for edges e, f is {a e ,aj} = 2 ^ cusps w({aj}, {&.,•}), for 
{oj}, {bj} the weights for the shears as sums of lengths of ideal geodesies. 
The matter is to evaluate the sum ([7]) for uj for the possible configurations. 
We first consider the case of the quadrilateral for the side e embedded in the 

2 We have reversed Penner's original sign convention given that his bivector has sides 
enumerated in a clockwise order, while Fock's bivector has sides enumerated in a counter- 
clockwise order. 
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surface and then describe necessary modifications for sides of the quadrilat- 
eral coinciding. The quadrilateral with weights for the edge e is given in 
Figure 




Figure 7: The quadrilateral for a triangulation edge e following formula 
([5]). The quadrilateral sides are labeled by lower case letters and vertices 
are labeled by Roman numerals. The edge weights 0, ±1 refer to expressing 
the e shear coordinate as a balanced sum of edge lengths. The numbers in 
square brackets are the sums Aj + Aj-\. 

Referring to formula ([7]), the first calculation is for the partial sums Aj 
of edge weights. At a vertex, edges are enumerated for summation in the 
counterclockwise order with the first edge being the clockwise most edge. 
Normalize the partial sums to be zero for the not listed edges preceding the 
first edge. The partial sums by vertex and in counterclockwise order are 
given in Table [TJ The second calculation is for the sums Aj + Aj-i of partial 
sums about vertices. The sums are given in Figure [7J by the numbers in 
square brackets; again sums vanish for edges not listed. Now we are ready 
to consider the configuration of the quadrilateral for the edge / and the 
sum of weights \{Aj + Aj—i)bj. The weights for / are again 0, ±1 as in 
Figure [7J The edges e and / are necessarily distinct. First consider that 
/ coincides with a boundary edge of the e quadrilateral. In this case the 
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Vertex 


Ai 


A 2 


^3 


I 


-1 


-1 





II 


1 







III 


-1 


-1 





IV 


1 








Table 1: Partial weight sums in counterclockwise order about vertices. 



diagonal edge weight for / is multiplied by the [±1] boundary edge weights 
for e and added to the ±1 boundary edge / weight times 1/2 the sum of 
the [2] and [2] diagonal weights for e. The result is ±2 with the positive 
sign if / is counterclockwise from e. Now consider the case that the e and 
/ quadrilaterals are either disjoint or intersect along a boundary edge. In 
the case of intersection along a boundary, the vanishing sum [1] + [—1] of 
e boundary weights gives a vanishing overall contribution. This completes 
the calculation if the quadrilateral of e is embedded. 

In general a pair of sides of the quadrilateral of e could coincide; we do 
not consider the special cases (g,n) = (0,3) or (1,1) where two side pairs 
coincide. A pair of adjacent sides could coincide by a 3/4 rotation about 
the common vertex or opposite sides could coincide by a translation. When 
sides coincide the contribution to lo is found by adding the contributions 
from each of the relative configurations for the quadrilateral of /. The result 
will be 0, ±4 according to adjacent or opposite sides coinciding and the e, / 
orientation. As already noted, we are using the adjustment jj Wlp07, §5] to 



our formulas 2luwp( ,t*) = d£* in place of ujwp{ ,t*) = d£* systematically 
used by Penner and Fock. The consequence is that our shear pairing is 
fourfold the Fock and Penner calculations. With this information, the shear 
pairing evaluations correspond and the proof is complete. 

□ 



9 The norm of a length gradient for a collar cross- 
ing geodesic 

We continue to consider compact surfaces with crossing geodesies a and 
a reflection symmetry, see Figure [3j We consider surfaces R t obtained by 
doubling a surface with cusps with ideal geodesies a, and opening cusps to 
obtain short length core geodesies 7. We are interested in the products of the 
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gradients grad£ a and grad£ 7 . Theorem [3] and Lemma [20] can be combined 
to provide expansions for the pairings 



Considerations of Chatauby convergence and sums of the differential f2 2 
from Section 2 suggest the heuristic expansion grad^ a = c a (£ 7 )grad^ 7 + 
ipit^ids 2 )- 1 with ip(ly) G Q(Re) converging to ip{0) G Q(RUR). A simple 
argument provides that ip(0) is orthogonal to the limit of grad l~. The above 
pairing formulas and heuristic then suggest an expansion 



The divergence of the pairing corresponds to the geometry. The limit of d£ a 
is formally the differential of length of an ideal geodesic and is a holomorphic 
quadratic differential with double poles at cusps. The limit is not an element 
of Q(R). Also the limiting infinitesimal deformation grad^ Q corresponds to 
opening cusps and has infinite WP norm. 

We would like to now use the gradient pairing formula, Theorem EJ 
to find the WP pairing for balanced sums of gradients of lengths of ideal 
geodesies. The above considerations show that a pairing formula involves 
canceling divergences in £ y . The divergences appear directly in evaluating 
the formula. The crossing geodesic a is orthogonal to the collar core 7. 
Arcs along 7 connect the intersection points with a. Each connecting arc 
provides a summand for the Theorem [3] evaluation. The connecting arcs 
along 7 occur in families; a family consists of a simple arc and the additional 
arcs obtained by adjoining complete circuits of 7. With l-y tending to zero 
and the summand i?(cosh dist) ~ 21og2/dist for small distance, there is an 
immediate divergence. We consider the sequence of lengths as a partition 
for a Riemann sum and find the £ 7 -asymptotics of the sum. 

The resulting formulas involve an elementary function, a reduced length 
for an ideal geodesic and a reduced connecting arcs sum formula. 

Definition 30. ForO <a<l, define the function A(a) = a(l — a)/(2sin7ra) 
with value given by continuity at the interval endpoints. For a crossing 
geodesic a on a compact surface R with reflection symmetry, the reduced 
length red(^ Q ) is the signed length of the segment connecting length 1 bound- 
aries of the complement of collars about core geodesies. For an ideal geodesic 




and 



(grad £ a , grad £ 7 ) = — (#a n 7) + 0(£). 



(grad £ a , grad £ a ) = — (#an 7 ) 2 + O(l). 
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a on a surface with cusps, the reduced length red(£ a ) is the signed length of 
the segment of a connecting the length 1 horocycles about the limiting cusps. 

The function A(a) is symmetric about a = 1/2 and satisfies 1/8 < A < 
1/27T. For a pair of points p, q on a circle, we write X(p, q) for the evaluation 
using the fractional part of the segment from p to q. For a hyperbolic surface 
without cone points the length 1 horocycles are embedded circles bounding 
disjoint cusp regions and red(^ Q ) is non negative. For surfaces with cone 
points, the reduced length can be negative. 

For crossing geodesies a, ft on a surface with reflection symmetry or ideal 
geodesies a, ft on a surface with cusps, we will write 



for the reduced sum over homotopy classes rel the closed sets a, ft of arcs 
connecting a to ft, that are not homotopic to arcs along a core 7 or along a 
horocycle. For the double of a surface with cusps, the symmetric homotopy 
classes are even with respect to the reflection; for this situation the sum is 
only over arcs with representatives on a chosen side of the surface. Each 
geodesic representative for the reduced sum intersects the thick subset of 
the surface and the reduced sum includes any intersection points of the ideal 
geodesies a and ft. We assume the main result Theorem [32] and illustrate the 
approach with the example of a single core geodesic. The general formula 
depends on the pattern of crossing geodesies. 

Example 31. Expansion of the WP gradient pairing for crossing geodesies 
a, ft and a single core geodesic 7. For the core intersections a n 7 = 
{a\, 02}, ft n 7 = {b\, 62} an d a given positive constant c then 

(grad^grad^,) = 




K 



2("> + red( ^) + 4 + 2 J- logA(a i) a i )) + K + 0(1^) 



(0i,aj) 



and for a ^ ft 



(grad 4,, grad^j) 
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We are ready to consider that pairings of balanced sums on a surface 
with cusps are the limits of pairings of balanced sums on approximating 
symmetric compact surfaces. The balanced sum condition will serve to can- 
cel the universal 16/£ 7 leading divergence terms. To compare formulas note 
that a surface with cusps represents half of a compact surface. It is also 
important that remainder terms as in the example tend to zero with . 

We state the main result. For a surface with cusps, the sum over core 
geodesic intersections is replaced with a double sum. First, a sum over cusps 
and second, a sum over ordered pairs of ideal geodesic segments limiting to 
a cusp. Ideal geodesies are orthogonal to horocycles. The fractional part 
of a horocycle defined by a pair of ideal geodesies is independent of the 
choice of horocycle. The geometric invariant A is evaluated by considering 
the intersections with any horocycle for the cusp. We present the formula 
for the case of a torsion- free cofinite group. 

Theorem 32. The ideal geodesic complex gradient pairing. For a surface 
R with cusps and balanced sums A = Yl a j^aj ,B = ^2 ^k^/3 k of ideal geodesic 
length functions, the WP pairing of gradients is 

(gradL(^),gradL(S)) = 

(2 2 
8 a . Ph -(ved(£ aj ) + 2) + - Yl ]og Xi&jJk) 

J> K cusps segments ctj,/3k 

limiting to the cusp 

n). 

ay to At J 

The first sum is over weights; the double sum is over ordered pairs of geodesic 
segments limiting to cusps. The final sum is over homotopy classes rel the 
closed sets ctj,/3k of arcs connecting ctj to (3k, arcs that are not homotopic 
into a cusp. For the homotopy class of an intersection ay n (3k, the function 
1Z is evaluated on cos 9, 9 the intersection angle. Otherwise, the function 
TZ is evaluated on the hyperbolic cosine of the length of the unique minimal 
connecting geodesic segment. Twist-length duality and J an isometry provide 
that 4(o'Ai ~b} = (grad^4, gradB). 

Proof. Begin the consideration with compact surfaces with reflection symme- 
tries and balanced sums of geodesic-length functions converging to a surface 
with cusps formally doubled across the cusps. The approach is to show that 
the connecting arcs sums of Theorem [3] converge to the sum for the limiting 
surface. The individual summands are considered in terms of the geometry 
of the biorthogonal connecting geodesic segments. 
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Begin by normalizing the uniformizations to ensure Chabauty conver- 
gence of the deck transformation groups T. For the geodesic a, let a be a 
chosen geodesic line lift and (A) the cyclic group stabilizer. A fundamental 
interval on a is chosen; each left (A) orbit in Ta and Tf3, f3 a lift of f3, has a 
unique biorthogonal geodesic connecting segment with one endpoint in the 
a fundamental interval. The considerations proceed in terms of the geom- 
etry of the second endpoint of the connecting segment. The finite number 
of terms corresponding to endpoints in a given compact set converge. The 
sums for families of connecting segments along the core geodesies provide 
universal divergences; the analysis is described in the next section. The 
remaining connecting segments have second endpoint outside a given com- 
pact set and the segments do not lie along core geodesies. The remaining 
segments necessarily intersect the lift of the thick subset. The remaining 
segments are treated according to whether the second endpoint lies in the 
lift of the thick or the thin subset. In the first case, the injectivity radius is 
bounded away from zero and the sum of such terms is uniformly bounded 
by applying the distant-sum method of [WlplO Chap. 8]. In the second 



case, the endpoint lies in the lift of a standard collar or cusp region. Hyper- 
bolic geometry is used to show that the full sum over the stabilizing cyclic 
hyperbolic or parabolic group is bounded simply by the distance of the fun- 
damental interval on a to the boundary of the region. The distant-sum and 
cyclic group bounds provide that the contributions from the complement of 
a large compact set is sufficiently small. The estimates for the various cases 
are combined to establish convergence of formulas. 

We consider the connecting segments along a given core geodesic. We 
outline the approach and give a detailed treatment in the next section. The 
sum for a family of connecting arcs in a given direction along a core geodesic 
has the form 

oo / 1 . , -I 

YS((a + n)i) for Sit) = cosh t ( log C °\ 
^-^ \ cosh t — 1 

n=0 v 

for t the core length and al, a > 0, the distance between core intersection 
points. The function S(t) has the initial expansion S(t) « 21og2/£ and for 
N approximately e > 0, we break up the sum 

oo 

J2S((a + n)e) = 

n=0 

N 2 1 N f 2\°° 

n=0 v ; n=0 v \ ' / / n =N+l 
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For the first sum, we use additivity of the logarithm to obtain an expression 
in terms of log 2jl and log T(a+1) for the gamma function. Stirling's formula 
is then applied. For the second sum, half of the first and last sum terms are 
separated, then the Trapezoid Rule is applied to approximate the sum by 
an integral and an error term. The Trapezoid Rule provides an improved 
approximation in i. The integral is calculated by an antiderivative. Finally 
the bound that S(t) is 0(e~ 2t ) for t > to > 0, provides that the third sum is 
exponentially small; the consequence is that for a > the original full sum 
has the expansion 

The overall expansion for connecting arcs in the forward and reverse direc- 
tions is obtained by combining the expansions for the values o and 1 — a. 
Identities for the gamma function are used to simplify the resulting formula 
and to obtain the function A. As already noted, the ^-divergence is in the 
leading term. The balanced sum condition provides for the overall canceling 
of divergences in evaluating the gradient product. The proof is complete. □ 

Corollary 33. For a balanced sum A = cij£ aj of ideal geodesic length 
functions and f3 a closed geodesic, the shear and twist derivative pairing is 

OAlp = -tpL(A) = ^Jdj cos6»p 

for the intersection angles measured from atj to f3. 

Example 34. A distance relation for the elliptic modular tessellation. 

The Dedekind tessellation is the tiling of the upper half plane for the 
action of PSL(2; Z). The light, respectively dark, triangle tiles form a single 
PSL(2;'Z) orbit. The reflection in the imaginary axis normalizes the group 
and interchanges the light and dark triangles. The tessellation vertices are 
fixed points of elements of the group action. There are two orbits for vertices. 
There are also two orbits for ideal lines. The first consists of the lines 
containing a single order-2 fixed point. The second consists of the lines 
sequentially containing an order-3, an order-2 and an order-3 fixed point. 
We refer to the types as 2-lines and 323-lines. We consider the lines with 
weights: w = +1 for 323-lines and w = —1 for 2-lines. The system of 
weighted lines is PSL(2; Z) invariant. 

The formula of Theorem [32] provides a relation for the distances between 
lines for the Dedekind tessellation. For any choice a of a 323-line and a of 
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Figure 8: The Dedekind tessellation. Graphic created by and used with 
permission from Gerard Westendorp. 



a 2-line, we have 



E 



w{rj)R(d{a,rjj) 



E 



w(r])R(d(a, rj)) 



log 



226 



ultraparallels to a 



ultraparallels to a 



for R(d) = it log ((u + l)/(u — 1)) — 2 and u = coshd. Ultraparallels are the 
tessellation lines at positive distance. Lines at zero distance are asymptotic. 

We find the relation as an exercise in evaluating the formula of Theorem 
[32l We begin with the geometry of the tiling quotient. We work with the 
thrice-punctured sphere uniformized by the projectivized index 6 subgroup 
PT(2) C PSL(2; Z) of matrices congruent to the identity modulo 2. A 
fundamental domain for the torsion- free group PT(2) is given by the twelve 
light and dark triangles adjacent to a given largest height non vertical 323- 
line. The PT(2) quotient is a tri-corner pillow with three 323-lines, labeled 
a,b,c and three 2-lines, labeled a, /3, 7. The 2-lines separate the quotient 
into two ideal triangles. A 323-line enters a single cusp of the quotient, 
while a 2-line connects two distinct cusps. We evaluate the pairing product 
for the weighted balanced sum a = a + b + c — a — f3 — 7. The sum is PT(2) 
invariant, thus grade E Q(PT(2)) by Theorem 12 li The space of PT(2) 
quadratic differentials is zero dimensional. The self pairing of grad a is zero. 

We determine the contributions for terms on the right hand side of the 
Theorem [32] formula. The evaluation corresponds to the formal expansion 
of the product (a + b + c — a — (3 — j) 2 . The pairing is real and the initial 
factor tt/2 can be moved to the left hand side. We begin with the reduced 
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length contribution. The PT(2) cusps have width 2; the length 1 horocycle 
at infinity has height 2. For a vertical 2-line, half of the reduced length 
segment connects the height two horocycle to the order-2 fixed point at 
height 1. A 2-line has reduced length 2 log 2. For a vertical 323-line, half of 
the reduced length segment connects the height two horocycle to the order-2 
fixed point at height 1/2. A 323- line has reduced length 4 log 2. The reduced 
length contributing terms of the product are a 2 + b 2 + c 2 + a 2 + /3 2 + 7 2 . 
The total first term reduced length contribution is 18 log 2 + 12. We next 
consider the log A contributions, which measure the geometry of the ideal 
geodesies limiting to cusps. There are two reflections stabilizing each cusp. 
The reflections stabilize the geodesies and provide that the intersections of 
the ideal geodesies with a horocycle are equally spaced and alternate by 
weights. The log A contributing terms of the product are 

a 2 + b 2 + c 2 + a 2 + f3 2 + 7 2 - 2a/3 - 2a 7 

- 2ba - 26 7 - 2ca - 2c/3 + 2a/3 + 2a 7 + 2/3 7 . 

By PSL(2; Z) symmetry, the evaluation is the same as for 3a 2 + 3a 2 -12a/3+ 
6a/3. The a 2 contribution is 2 log(A(0)A(l/2)) given the two segments at a 
cusp; the a 2 contribution is 21ogA(0) given the two limiting cusps; the a/3 
contribution is 2 log A (1/4) given the symmetry of A and the a/3 contribution 
is logA(l/2). The evaluations are A(0) = l/(2vr), A(l/4) = 3a/2/32 and 
A(l/2) = 1/8. The total log A contribution is 

1 1 3^/2 1 

61 ° g 16^ + 61 ° g 2^ + - 24l0g ^2" + 61 ° g 8- 

We next consider the contribution from ideal geodesies intersecting. The 
intersection product contributing terms are 2ab+2ac-\-2bc — 2aa — 26/3 — 2c 7 . 
The geodesic intersections ab, ac and be are twofold. From the formula the 
total intersection contribution is 

2 • 3 • i?(cos -) + 2 • 3 • R(cos — ) - 2 • 3 • ii(cos -) = 6 log 3 - 12 

as follows. The leading 2-factors are from the formal expansion of a 2 . The 
3-factors are from the symmetry of the triples a,b,c and q,/3, 7 . The first 
and second terms correspond to the fact that distinct 323-lines intersect 
twice. The ^-evaluations i?(cos7r/3) = (log3)/2 — 2 and i?(cos7r/2) = —2 
are elementary. The final contribution of the right hand side of the overall 
formula is the sum for the nontrivial connecting geodesies. We start with the 
formal expansion a 2 = aa + ba + ca — aa — f3a — , ya. By PSL(2\7L) symmetry 



47 



the evaluation is the same as for 3aa — 3aa. Connecting geodesies are 
enumerated by lifting to the universal cover. Given lifts a and a, the desired 
sums are obtained. The overall relation now follows. We note that the lines 
asymptotic to a and a correspond to the limits of lines with connecting 
segments along core geodesies; the log A terms account for the combined 
contribution of the asymptotic lines. 

10 The geodesic circuit sum 

We consider the contribution to the Theorem [3] sum corresponding to con- 
necting geodesies given by circuits about a fixed closed geodesic. Such a 
circuit sum enters when the geodesies a and f3 are orthogonal to a common 
closed geodesic. The summands are evaluations of the function 



The consideration is for the length parameter I expansion of the infinite 
sum of circuits. The application to Theorem 1321 requires an expansion with 
remainder term tending to zero for small I. Simple analysis gives that the 
expansion begins with terms divergent in t. We provide the expansion. 

Theorem 35. For a and e positive, the circuit sum has the expansion 



for the gamma function T(z). 

Corollary 36. For e positive, the circuit sum for < a < 1 has the expan- 
sion 




^5((a + n)£) = - +log 



r(a + l) 2 ^ 2a " 1 



+ 2a - 1 + Oil 1 ' 6 ) 



n=0 



S((a + n)£) = - + 2 log 



T(a + 1)T(2 - a) 
2^ 



+ oil 1 - 6 ) 



n=— oo 




and for a = 1 has the expansion 



2 f 

J2s(n£) = - + log + l + 0(£ 1 -*). 



n=l 
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Proof of Corollary. Since S(t) is an even function the first sum can be 
rewritten as YlnLo S(( a + n K) + &(0 ~ a + n K) an< ^ * ne theorem is applied. 
The gamma function identities T(z + 1) = zT(z) and T(l — z)T(z) sin7rz = ir 
are applied to obtain the desired expression. Finally the case a = 1 is a 
direct application of the theorem. □ 
Proof of Theorem. We begin with properties of the summand S(t). The 
summand has the small-i expansion S(t) = 21og2/t - 2 + 0(t 2 logt) and 
the large-t expansion S(t) = 0(e~ 2t ). We also consider the function 

F(t) = S(t) - 2 log - t 

and write 

^, s , , J, coshi+l\ , t 2 (cosht + l) 

F(t) = (cosht-1 log + log -7 — : - - 2. 

w v \ cosht-iy & 4(cosht-l) 

We note that for small-t, since cosht — 1 is 0(t 2 ) and £ 2 /(coshi — 1) is 
analytic it follows that F(t) has second derivative bounded by — logi for 
small-i. 

We are ready to begin the overall considerations and write the sum in the 
form of Riemann sums, adding in and subtracting out a 2 log 2/t contribution 



2 



n=0 n=0 v ; 

1 N 2 
+ -^4% + n W -21og— -) 



e^L V VV ' ' G (a + n) 

n=0 v ' 

CO 

+ 7 E ^((a + n)t) 

n=N+l 

= i + n + in. (9) 

We consider the right-hand sums in order. For the first sum we have 

N 2 2*1 

%^UTny = + + %^Jar+-nY 

n=0 v ' n=0 v ' 

The right hand sum is -2 log n^ =0 (a + n) = -2 log Y(a + N + l)/r(a + 1). 
We apply Stirling's formula logr(z) = ilog27r/z + z(logz — 1) + 0(1/ z) 
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to find that 

/ = 2(JV + l)log 7 + 21o g r(a + l) - 2(a + N + ]-)log{a + N + 1) 

+ 2(a + iV + l) - log 2tt + C^iV" 1 ) 

and noting that log (a + N + 1) = log (a + N) + I /{a + N) + 0(N~ 2 ) gives 
the desired final expansion 

/ = 2(iV + l)log 7 + 21o g r(a + l) - 2(a + N) log(a + iV) 

- log(a + iV + 1) + 2(a + JV) - log 2vr + 0(iV _1 ). (10) 

For the second sum of (jUJ) we use the Trapezoid Rule approximation 
for an integral. The approximation involves weights 1/2 for the first and 
last sum terms. The error bound is in terms of the second derivative of 
F(t) on the interval [a£, (a + N)£] and the square of the partition size. The 
approximation gives the expansion 

i r(a+N)l i 

11= j J F(t)dt + -{F(a£) + F((a + N)£)) 

+ 0(£\[a£,(a + N)£]\max\F"\). 

We set (a + N) = £~ e and consider terms in order from right to left. Given 
the small-f logarithmic bound for F" the remainder is bounded as 0{£ l ~ 2)L ). 
Given the large-i exponential decay S(t) and the small-t expansion of S(t) 
then 

F((a + N)£) = -2 log * + 0{e~ r ') and F(a£) = -2 + 0{£ 2 ~% 

The next step is to include the contribution of sum III. The sum is re- 
placed with the corresponding integral. Since the integrand is exponentially 
decreasing on the interval, the replacement remainder is exponentially small. 
The considerations combine to give the expansion 

9 p(a+N)£ 9 1 f 00 

11 + 111 = — log dt + - S(t)dt 

1 J at 1 1 J at 

- 1 - 1 ° g (^K + 0(fI " , »- 
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The first integrand has antiderivative ilog2/i + t. The second integrand 
S(t) has antiderivative 



. . , / . cosh t + 1 
smhi log — — — - 
\ cosn t — 1 

which has the large-t expansion 2 + 0(e -2 *). We evaluate the integrals to 
find the contribution 

II + /// = -2iVlog2 + 2(a + N)log((a + N)£) - 2(a + N) 

- 2alo S at + *,+ \ - 2alog * - 1 - log,;^ + 0(1^). 



The next step is to combine with expansion (llOj) and note again that (a + 
N)£ = £~ e to find the desired final expansion 

I + II + III = j + log lQ + 2a - 1 + O^ 1 -). 
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